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This paper deals with the large-time analysis of a PDE system modelling contact with 
adhesion, in the case when thermal effects are taken into account. The phenomenon of adhe- 
sive contact is described in terms of phase transitions for a surface damage model proposed 
by M. Fremond. Thermal effects are governed by entropy balance laws. The resulting 
system is highly nonlinear, mainly due to the presence of internal constraints on the phy- 
f^T) ' sical variables and the coupling of equations written in a domain and on a contact surface. 

ON , We prove existence of solutions on the whole time interval (0, +00) by a double approxima- 

tion procedure. Hence, we are able to show that solution trajectories admit cluster points 
which fulfil the stationary problem associated with the evolutionary system, and that in the 
ON ' large-time limit dissipation vanishes. 

o 
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■ - - ■ 1 Introduction 



This paper is concerned with the large-time analysis of a PDE system describing adhesive contact 
between a thermo-viscoelastic body and a rigid support. The model has been recently introduced 
and global-in-time existence results have been proved on finite-time intervals both in isothermal 
cases (see [2] in the case of an irreversible damage evolution on the contact surface, and [3] for 
the reversible case), and for PDE systems including thermal effects (see [1]). The modelling 
approach for contact with adhesion which we apply refers to a damage theory described by 
phase transitions, and it is due to M. Fremond (see |13j). The idea consists in describing the 
adhesion between viscoelastic bodies in terms of a surface damage theory, in which the damage 
parameter is related to the active bonds which are responsible for the adhesion between the 
bodies. Hence, the equations of the evolutionary system are recovered from thermomechanical 
laws and they are written in the domain of the viscoelastic body and on the contact surface. 
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It turns out to be interesting, both from a theoretical point of view and in view of applications, 
to investigate how the thermomechanical system (i.e., the body and the rigid support it is in 
contact with) behaves for large times. More precisely, we shall investigate if the trajectories 
of the solutions to the resulting PDE system present some cluster point, in the limit as time 
goes to +00. Then, we shall look for a relation between these limit states and the stationary 
system associated with our evolution problem. In particular, we aim to prove that, in the limit, 
solution trajectories reach a thermomechanical equilibrium state in which dissipation vanishes. 
This kind of large-time analysis was performed in [3] for the reversible model in the isothermal 
case. 

Before introducing the long-time behaviour analysis of the problem, we shall briefly recall the 
model and make some comments on the existence of solutions on the whole time interval (0, +00). 
Moreover, we shall point out that this paper also presents a novelty in the formulation of the 
model itself, as we generalize the convex potential usually ensuring internal constraints on the 
damage parameter. 

The model 

We mainly refer to the recent contribution [3], in which (a slightly different version of) the ther- 
momechanical model has been introduced. The state variables, in terms of which the equilibrium 
of the system is established, are defined in the domain ft C ffi 3 (where the body is located), 
and on the contact surface T c . Namely, we shall take f! to be a sufficiently smooth bounded 
domain in M 3 , with boundary dfl = f 1 U T2 U F c . The sets Tj are open subsets in the relative 
topology of 90, with smooth boundary and disjoint one from each other. In particular, T c is the 
contact surface. We suppose that T c and T\ have strictly positive measures and, for the sake of 
simplicity, we identify T c with a subset of M 2 , i.e., we shall treat T c as a flat surface. 

The state variables we shall consider in Q are the absolute temperature $ of the body and 
the macroscopic deformations, given in terms of the linearized strain tensor e(u) (u represents 
the vector of small displacements). On the contact surface, we introduce the surface absolute 
temperature (the reader may think of the temperature of the adhesive glue) i? s and a damage 
parameter X, related to the active bonds in the glue ensuring adhesion. For the moment, we do 
not require any constraints on the values assumed by X. Taking into account local interactions 
(in the glue and between the glue and the body) we include the gradient VX and the displacement 
trace U| r among the state variables on the contact surface. The free energy in £1 is written as 
follows 



where K is the elasticity tensor and the coefficient *& multiplying tre(u) accounts for the thermal 
expansion energy. 

Remark 1.1. Notice that here we have taken the term $(1 — ln(i9)) for the purely thermal 
contribution in the free energy ^/q (and, similarly, for the free energy ^r c below), while in 
[1] we have considered a more general concave function. The particular choice in this paper 
is very frequent in the literature, as it has some analytical and modelling advantages. From 
the latter viewpoint, the presence of the logarithm in yields an internal constraint on 

the temperature: indeed, the domain of is given for 1? > 0, which is in agreement with 
thermodynamical consistency. On the analytical level, this form of the thermal contribution 
shall allow us to simplify the procedure exploited in [3J to prove existence of solutions (see 
Remark S3]). 
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Next, we specify the free energy in T c , which presents some novelty with respect to the model 
introduced in [1] (cf. also [2] and [3]). In fact, we shall consider 

*r c = 0(1 - ln(^)) + A(X)(0. - V eq ) + (3(X) + a(X) + ||VX| 2 + V|u, rc j 2 + /_( U | rc • n), (1.2) 

where $ eq is a critical temperature, and (3 is a proper, convex, and lower-semicontinuous function. 
The indicator function 2_ forces the scalar product U| r • n to be non-positive, as it is defined on 
K by I-(y) = if y < and I—(y) = +oo for y > 0. This renders the impenetrability condition 
between the body and the support. In the same way, the term (3 may yield a constraint on 
the values assumed by X. For example, a proper choice of (5 may enforce positivity of X (see 
[21 El H]). In particular, this occurs when, classically, = I[oi], forcing X € [0,1]. Then, the 
coefficient of |u| r | 2 remains non-negative, in accord with physical consistency. However, in the 

present paper we shall allow the potential (3 to be more general and we do not impose any a 
priori restriction on its domain. Hence, to ensure physical consistency, we introduce a constraint 
on the deformation coefficient, which is fixed to be X + (using the notation r + = max(r, 0) for 
every r£l). Finally, the function A is related to the latent heat, while a takes into account 
possibly non-convex contributions in the free energy. In particular, we include in a cohesive 
effects in the glue, which are represented by a non-increasing function in X (a simple choice is 
cr(X) = w(l — X), with a positive parameter w). 

Then, the evolution of the system is governed by two convex potentials (non-negative and 
assuming their minimum if there is no dissipation) , namely the pseudo-potentials of dissipation 
written in and in r c . We have (here K v is a viscosity matrix) 

<5>n = \\V$\ 2 + \e{vit)K v E{u t ), (1.3) 

and 

*r c = ^|W S | 2 + \\X t \ 2 + \k{X){\ c - d s f. (1.4) 

Notice that the dissipation in depends on e(uj) and on Vi?, while the dissipation in T c depends 
on Vi? s , Xt, and on ($i r — i? s ). The function k, which accounts for the heat exchange between 
the body and the adhesive material, shall be taken non-negative and smooth enough. Actually, 
to characterize the large-time behaviour of the system, we need to assume that k is bounded 
from below by some positive constant (see Remark II. 3p . 

The PDE system 

Proceeding as in we refer to thermomechanical laws and, after specifying the constitutive 
equations in terms of the above potentials, we arrive at the following PDE system (T is a fixed 
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final time) 

dt(hx(0)) - div(ut) -A# = h in Q x (0,T), (1.5) 

(0 in(^\r c )x(0,T), 

\-fc(X)(0-0 8 ) inr c x(0,T), 

$(ln(0,))-3t(A(X))-A0 s = fc(X)(0-i?,) inr c x(0,T), (1.7) 

d n $ s = mdT c x(0,T), (1.8) 

- div (ife(u) + ^e(ut) + 01) = f in x (0, T), (1.9) 

u = inrix(0,T), (ife(u) + K v e(u t ) + 01)n = g inr 2 x(0,T), (1.10) 

(Ke(u) + K v e(u t ) + 01)n + X+u + <9/_(u • n)n 3 inr c x(0,T), (1.11) 

X t -AX + f3(X) + o-'{X)-\'(X)$ eq + R(X)^\u\ 2 B-\'(X)$ s mT c x(0,T), (1.12) 

a ns x = o mar c x(o,r) (1.13) 

where /i is an external entropy source, f a volume force, and g a traction. Moreover, (3 = d(3 
and H = dp, where p(r) = r + for all r£l. Hence, the Heaviside maximal monotone operator 
H : R -> 2 R is defined by H(r) = if r < 0, H(0) = [0, 1], and H(r) = 1 if r > 0. We warn that, 
here and in what follows, we shall omit for simplicity the index v\ r to denote the trace on T c 
of a function v, defined in fl. 

Remark 1.2. Let us comment on the above equations, while referring the to [I] for their rigorous 
derivation. First of all, we point out that (|1.5p and (|1.7|) are entropy equations. The possibility 
of describing thermal effects in phase transitions by the use of an entropy equation, in place of 
the more standard energy balance, has only recently been introduced. In particular, let us point 
out that the entropy in O is defined as ln(0) — div(u) and in T c as ln(0 s ) — X(X). Using entropy 
equations brings to some advantages both for the analytical treatment and the modelling of 
the phenomenon. In particular, from (ll.5p and (II. 7h one directly recovers the positivity of the 
temperature, which is necessary for thermodynamical consistency, avoiding the application of 
any maximum principle argument. We do not enter the details of this theory and refer, among 
the others, to the papers [7] and [8]. Then, (|1.9p is derived from the momentum balance, in which 
accelerations are not taken into account. Equation (|1.12p is recovered as a balance equation for 
micro-movements related to the evolution of the phase parameter (see [13] for the theory of the 
generalized principle of virtual power including micro-movements and micro-forces responsible 
for the phase transition). 

Remark 1.3. We emphasize that the structure of (|1.12p is more complicate than the analogous 
equation in the model studied in [1]. Indeed, the maximal monotone operator /3 is more general 
than the one considered in [3], for we do not impose any restriction on its domain. Moreover, 
the presence of the operator H in (I1.12|) introduces a new nonlinearity in the equation. From 
the physical viewpoint, since H(0) = [0, 1], a residual influence of macroscopic displacements on 
the mechanical behaviour of the glue may persist also when the glue is damaged, e.g. if X = 
(X representing here the proportion of active bonds). This corresponds to assuming that a local 
interaction between the body and the support is preserved even when the bonds in the glue are 
completely damaged. Notice that this is reasonable, if one takes into account distance forces. 
An analogous argument justifies the assumption that k in (|1.7p is bounded from below by some 
positive constant, see ()2.H5p later on. This ensures that the thermal local interaction between 
the body and the support is conserved even when the adhesion is not active. 
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Our first main result (see Theorem 12. II later on) states that for every T > the Cauchy problem 
for system (|1.5p - (|1.13|) admits at least one solution. In this way, we parallel the global existence 
result of [3]: therein, as we mentioned before, we considered a slightly different free energy on the 
contact surface, which resulted in an equation governing the evolution of the parameter X simpler 
than (|1.12l) . Nonetheless, the proof of Theorem 12.11 (which shall be developed in Section 0]), 
closely follows the argument developed in [3]. It hinges upon a double approximation procedure 
(depending on two approximating parameters), and a subsequent passage to the limit argument 
with respect to the mentioned parameters. One of them is used to regularize the nonlinearities 
in the equations by means of Yosida approximations. Furthermore, some viscosity terms in i? 
and $ s are added in (|1.5p and (|1.7p . depending on the second parameter. The local existence 
of a solution for the approximate system (supplemented with suitable regularized initial data 
for ■& and due to the presence of viscosity), is obtained with the Schauder theorem, while 
uniqueness follows by contracting arguments. Hence, we conclude the existence of global-in-time 
solutions by proving suitable a priori estimates (which in fact directly hold in the time interval 
(0, +oo), as they do not depend on the final time horizon T), independent of the approximating 
parameters. The very same estimates allow us to pass to the limit in the approximate problem, 
firstly as the viscosity parameter and secondly as the parameter of the Yosida regularizations 
vanish. Finally, we point out that, due to the strongly nonlinear character of system (|1.5p - (|l,13l) . 
we do not expect uniqueness of solutions for the related Cauchy problem. 

Large-time analysis 

As previously mentioned, the ultimate aim of this paper is investigating the large-time behaviour 
of system fjl.5p — fjl.l3j) (supplemented with suitable initial conditions). More precisely, we are 
interested in finding cluster points of solution trajectories and characterizing a sort of ther- 
momechanical equilibrium of the system in the limit, in which there is no dissipation. This 
corresponds to proving that solution trajectories converge to solutions of the stationary problem 
associated with our system, in which dissipation is zero. 

Now, some results in this direction have been obtained in the literature concerning the long- 
time behaviour of phase-field systems with non-convex potentials (see, for example, [111 [121 El 
[To] ll7j). Typically, these results apply to binary systems (i.e., macroscopic deformations are not 
included), see among the others [8] dealing with a singular entropy equation. 

The main difficulties related to our analysis are due to the singular character of the entropy 
equations, to the presence of general multivalued operators on the state variables, and to the 
nonlinear coupling between the equations written in the domain 0, and the ones set in r c . 

Remark 1.4. On the other hand, the analysis of the large-time behaviour in terms of the global 
attractor for the dynamical system generated by (|1.5|) - (|1.13|) might also be addressed. Indeed, 
the existence of the global attractor would signify that the system dissipation is controlled in the 
evolution. However, for the moment being, proving the existence of the attractor seems out of 
our reach. In fact, the strongly nonlinear character of the equations essentially prevents us from 
obtaining those estimates on the solutions which would guarantee the existence of a compact 
and absorbing set, for the dynamical system, in the phase space dictated by the choice of the 
initial data. 

Prior to addressing the large-time analysis of system (|1.5p - (|1.13p . we specify that we consider 
a quadruple (#,0 a ,u,X) to be a solution of (fT3I)- (fIT31) in (0,+oc), if (0,tf a ,u,X) fulfils ([13]) - 
(I1.13P in the finite-time interval (0, T), for every T > 0. Hence, to perform the asymptotic 
analysis on the solutions of (|1.5|) - (|1.13|) as time goes to +co, we shall rely on some further 
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estimates improving the solution regularity of the existence Theorem 12.21 Only in this enhanced 
setting, shall we obtain (see Proposition I2.5P the bounds on the solutions (in suitable functional 
spaces, on the whole half-line (0, +00)), necessary to prove that, for every solution trajectory, 
its cj-limit set (i.e., the set of its cluster points) is non-empty. These further estimates shall 
be first formally derived in Section [3l and then made rigorous in Section 0] by performing all 
the related calculations on the approximate system used for proving Theorem 12.11 That is 
why, our asymptotic analysis solely applies to the solutions of (|1.5jl - (|1.13p originating from the 
abovementioned double approximation procedure. Once proven that the u;-limit is non-empty, 
we shall show that its elements solve the stationary system associated with the evolutionary 
problem fjl.5j) — (|1.13j) (see Theorem 12.2 j) . As a by-product of this procedure, we shall see that in 
the limit as t — > 00 the dissipation vanishes (cf., in particular, Remark 12. 7D . 

Plan of the paper. In Section [2] we enlist all of the assumptions on the problem data and 
state our results. A (partially formal) proof of our Theorem 12.21 on the long-time behaviour 
of the PDE system (|1.5p - (|1.13|> is developed in Section [3] and rigorously justified in Section [H 
which also contains the proof of the global existence Theorem 12.11 

2 Main results 
2.1 Preliminaries 

Notation 2.1. Throughout the paper, given a Banach space X, we shall denote by x'i'i ')x the 
duality pairing between X' and X itself, and by || • ||x both the norm in X and in any power of 
X; C^([0, T];X) shall be the space of the weakly continuous AT-valued functions on [0,T]. 

Henceforth, we shall suppose that Q is a bounded smooth set of M 3 , such that T c is a smooth 
bounded domain of IR 2 , and use the notation 

H := L 2 (ty, V:=H\n), and 
W := {v € V 3 : v = a.e. on Ti} , 

the latter space endowed with the norm induced by V. We shall denote by 31 the standard Riesz 
operator 

%:V->V given by v , (Jl(u), v) v := uv + / VuVv for all u, v £ V, (2.1) 

Jo, Jn 

and by 0lr c the analogously defined Riesz operator, mapping if 1 (T c ) into (i? 1 (r c )) / . We shall 
extensively use that 

V C L P (T C ) with a continuous (compact) embedding for l<p<4 (1 < p < 4, resp.), (2.2) 

.ff 1 (r c ) C L P {T C ) with a compact embedding for 1 < p < 00. (2-3) 

For notational simplicity, we shall write J r uv uv, respectively) for the duality pairing 

(2f-Va(r c ))s(u,v) (H i/ a(ro))8 between (tf-^r^s an d (tf 1 / 2 ^)) 3 (between (H~ 1 / 2 {T 2 )) 3 and 

(i? 1 / 2 (r 2 )) 3 , resp.). Finally, given a subset C R N , N = 1,2,3, we shall denote by |0| 
its Lebesgue measure and, for a given v G V', the symbol m(v) shall signify its mean value 

i/\n\ v ,(v,i)v 
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Variational formulation of the elasticity equation. We introduce the standard bilinear 
forms which allow to give a variational formulation of (the boundary value problem for) equa- 
tion (|1.9p . As usual in elasticity theory, we may assume that the material is isotropic and hence 
suppose that the rigidity matrix K in (|1.9|) - (|1.11|) can be represented as 

Ke(u) = Atr (e(u)) 1 + 2pe(u), 

where A, /i > are the so-called Lame constants and 1 is the identity matrix. Also, for the sake 
of simplicity but without loss of generality, we set K v = 1 in (|1.9|) - (|1.11|) . Therefore, (|1.9|) may 
be formulated by means of the following bilinear symmetric forms a, b : W x W — > R, defined by 

a(u, v) := A / div(u) div(v) + 1\i \^ \ £ ij( u ) £ ij(' v ) f° r au u, v € W, 

Jn ij=1 Jn 

3 

i»(u, v) = 2_, / e ij( u ) £ ij( v ) f° r an u, v e W. 

iJ=i Ja 

Note that the forms a(-, •) and b(-, •) are continuous and, since T\ has positive measure, by Korn's 
inequality they are W-elliptic as well, so that 

3C a , K a >0: a(u,u) > C ||u||w |a(u, v)| < # a ||u|| w ||v|| w for all u, v G W, (2.4) 
3C b , K b >0: 6(u,u) > C 6 ||u||^ |6(u,v)| < J^||u|| w ||v||w for all u, v G W. (2.5) 

2.2 A global existence result 

Statement of the assumptions. In equation (jl . 121) we consider 

a maximal monotone operator (3 : R — > 2 K , (2. HI) 



and denote by (5 : D(/3) — > (— oo,+oo] a proper, l.s.c. and convex function such that (3 = d(3. 
Instead, of dealing with the pointwise operator dl- : R — > 2 K in (jl.lip . we shall work with a 
suitable generalization, defined in the duality relation between i? -1//2 (r c ) 3 and -ff 1 / 2 (r c ) 3 . To 
this aim, we introduce 

a : (i7 1 ^ 2 (r c )) 3 — > [0, +oo] a proper, convex and l.s.c. functional, 

^ ^ (2.H2) 

with 8(0) = = mina, 

and set 

a:=da: {H l l 2 {T c )f ^ 2^ H ' 1/2 ^ . 

Remark 2.2. In order to render the impenetrability constraint mentioned in the Introduction 
by means of the operator a, we may proceed as follows. We let j(u) = /_(u • n) and associate 
with j the following functional 

«(v) = / j(v) if v G (tf 1 / 2 ^)) 3 and j(v) G L 1 ^), (2.6) 

S(v) = +oo otherwise. (2-7) 

Since a is a proper, convex and lower semicontinuous functional on (i? 1 / 2 (r c )) 3 , its subdiffer- 
ential (cf. pQ Cap. II, p. 52]) 

a := da : (tf 1 / 2 ^)) 3 -> 2^ 1/2 ^ (2.8) 
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is a maximal monotone operator. Notice that, in this case, (|2.6[) implies that, if rj G ce(v), then 
v belongs to the domain of j and thus fulfils v • n < 0, which corresponds to the impenetrability 
condition. 

We assume that the nonlinearities a and A comply with 

ueC^R), (2.H3) 

AeC^R), (2.H4) 

(and denote by L a and L\ the Lipschitz constants of the functions a' : R — > R and A' : R — > R, 
respectively), and that (cf. Remark 1 1 .3D 

k : R — ► (0, +oo) is Lipschitz continuous , with Lipschitz constant Lk, and 

3c fc >0VxGR: k(x) > c k . (2.H5) 

Remark 2.3. We point out that ([27H3]) . ([2~H4]) . and ([2Jl5]) respectively entail that 



3C a > Vx G 
3C A >0Vxe 
aC fe > Vx € 



k(x)| <C CT (x 2 + l), (2.9a) 
|A'(x)| <C A (|x[ + l), (2.9b) 
\k(x)\ < C k (\x\ + 1) . (2.9c) 



As far as the problem data are concerned, we suppose 

h G L 2 (0, T; V) n 1^(0, T; fi) , (2.H6) 

fGL 2 (0,T;W), (2.H7) 

geL^O.T;^ 1 / 2 ^)) 3 ). (2.H8) 
It follows from (l2~H7l) - (l2~H8l) that, defining F : (0, T) -> W via 

w ,(F(t),v) w := w /(f(t),v) w + / g(t).v VvGty for a.e. t G (0, T), 

there holds 

F G L 2 (0, T; W') . (2.10) 
Finally, we require that the initial data fulfil 

i?o S LP(fl) , with p > -, and ln(tf ) G H , (2.11) 

5 

??° G L 9 "(r c ) , with g > 1, and ln(i9°) G L 2 (T C ) , (2.12) 

u G W and u G D(S) , (2.13) 

X G H\T C ) and 0(X O ) G L\T C ) . (2.14) 

Note that the first of (f2~PT|) and of (pA2|) respectively yield 

i?o^' ^ s eH\T c )'. (2.15) 
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Variational formulation and existence theorem. The variational formulation of the initial- 
boundary value problem for system (| 1 . 5 j) — (j 1 . 1 3 j) reads as follows. 

Problem (P). Under the standing assumptions (j2.Hlh ~ (j2.H8h . given a quadruple of initial 
data (i?o, uo, Xo) complying with (|2.11|) - (|2.14|) . find functions i? s , u, X, rj, £, (), with the 
regularity 



fulfilling the initial conditions 



and 



(2.16) 
(2.17) 



G L 2 (0, T; V) n L°°(0, T; L 1 (fi)) , 
ln(tf) G L°°(0, T; ff) n fl^O, T; V') , 
tf s GL 2 (0,T;i/ 1 (r c ))nL oo (0,r;L 1 (r c )), 
ln(0.) G L°°(0, T; L 2 (r c )) D fl^O, T; H 1 (T c )') , 

u G H\0,T;W), (2.18) 

T7GL 2 (0,T;(ff- 1 /2(r c ))3), (2.19) 

x g L 2 (o,T; J ff 2 (r c ))nL 00 (o,r;F 1 (r c ))n J ff 1 (o,r ; L 2 (r c )), (2.20) 

eGL 2 (0,T;L 2 (r c )), (2.21) 

CGL°°(r c x(0,T)), (2.22) 



!?(0) = tf a.e. in , (2.23) 

S (O) = 02 a.e. inr c , (2.24) 

u(0) = u a.e. infi, (2.25) 

X(0)=X o a.e. inr c , (2.26) 



(2.27) 



(2.28) 



(2.29) 



v ,(Qtln(i?),v) v - [ div(u t )v+ [ WVt> + ( k(X)($-# s )v 
Jn Jn Jr c 

= v ,(h,v) v V v G V a.e. in (0, T) , 
HH r c y{dtM'Vs),v) H i {Tc) - j d t \{X)v+f W,Vu 

= / fc(X)(0-0 s )u V^Gil 1 ^) a.e. in (0,T), 

6(u«,v) + a(u,v) + / tfdiv(v)+ / (X + u + r?)-v 

= w( F > v >w Vv G W a.e. in (0,T) , 
?7 G a(u) in (#~ 1/2 (r c )) 3 a.e. in (0,T), (2.30) 

X* - AX + £ + <7'(X) = -A'(X)t? s - ^C|u| 2 a.e. in T c x (0,T), (2.31) 

£ G (3(X) a.e. in T c x (0,T), (2.32) 
C G H(X) a.e. in T c x (0,T), (2.33) 
ns X = a.e. in dT c x (0, T) . (2.34) 

Notice that the contribution — A'(X)i? cq occurring in (jl . 121) has been incorporated into the term 
cr'(X) in (gj5ID . 
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Theorem 2.1. Assume (I2HT1) - rt27H8j) . 

1. Then, Problem (P) admits a global solution $ s , u, X, r], £, £) on i/ie interval (0, T). 

2. If, in addition, 

fG^TjW'), je^'fo.Tjfr 1 / 2 ^)) 3 ), /»€ w^co.TiV), (2.35) 

i/ien i/iere exists a solution u, X, 77, £, £) having for all 5 > t/ie further regularity 



€ L°°(o~, T; V) n ff 1 ^, T; L 12 / 7 (0)) , (2.36a) 

S GL°°(5,T;^ 1 (r c ))nii' 1 (^r;L 2 ^(r c )) /or aW p € (0,2), (2.36b) 

X G L°°(5,r;if 2 (r c )) n tf^T;^ 1 ^)) n^°°(^;L 2 (r c )) , (2.36c) 

^L°°(5,T;L\T C )), (2.36d) 

ueW^rjW). (2.36e) 



From (I2.36al )- (|2.36cp it follows in particular that for all 5 > 

# e C° W ([5,T};V) , $ S £C° W ([5,T];H\T C )), X G C° ([5, T}; H 2 (T C )) . (2.37) 



We refer to Remark 14.121 for some further comments concerning the above statement. 

Remark 2.4. The regularity of $o and $® required in (I2.1ip - (|2.12l) turns out to be necessary 
in the proof of our existence result Theorem 12.11 for (the Cauchy problem for) system (|1.5p 
(|1.13p . Indeed, since we are going to prove existence of solutions by passing to the limit in a 
viscosity approximation of equations (jl.5p and (jl.7p . we shall need to dispose of more regular 
approximate initial data. Our construction of such data (see Lemma I4.2|) apparently hinges 
upon the regularity fl2TTTD - (|2TT2] ) of tf and 0°. 

However, we are not able to recover for and S the regularity corresponding to assumptions 
(|2XTj) and (pT2l) . namely G C°([0, T\;L P (VL)) and S G C°([0, T]; L 9 ~(ft)), with p and q as in 
(I2.1ip - (|2.12|) . This is mainly due to the highly nonlinear character of PDE system and, in some 
sense, to the fact that the natural initial conditions for (|2.27p and (|2.28p are written for ln(0) 
and ln(0 s ). Nonetheless, notice that the regularity required for the initial data is preserved (see 
(I2.37P ) for t > 5 > 0, for every 5 > 0, in the more regular framework of (|2.35p . 

The proof of the above result is based on a double approximation procedure which we shall 
detail in Section [H The related passage to the limit relies on suitable a priori estimates on the 
approximate solutions, which we shall formally perform on the (un- approximated) system (|2,27|) 
(I2.34p . and directly on the time-interval (0, +oo), within the (formal) proof of Proposition 12.51 
Such estimates shall be rendered rigorous in Sec. 14.31 Finally, in Sec. 14.41 we shall conclude the 
proof of Theorem 12.11 

2.3 Results on the long-time behaviour of Problem (P) 

Within the scope of the present section, we shall say that 

a quadruple (0,0 s ,u,X), with the regularity (l2~TB - (l2~THI) and (|2~2"U1) . 
is a solution to Problem (P) if it fulfils fl2~2T)) - ([2~2"g]) and 
there exists a triple (»/,£, C) for which (I2TT9D . (I23TD . (I2T22D and (l2T29D - (l234l hold. 
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In view of the long-time analysis of the solutions to Problem (P), we shall hereafter suppose 
that 

Vi? > 3C R , > Vie dom(3) : &(x) +C R > Rx 2 . (2.H9) 

Notice that (|2.H9p is trivially fulfilled in the case dom(/3) is a bounded interval, whereas, if 
dom(/3) is unbounded, it is implied by a super-quadratic growth of (5 at infinity. We shall also 
require some summability on (0, +oo) for the problem data: 

f G L°°(0,+oo;W') and ft G L 1 (0, +oo; W') , (2.H10) 
ger^+oo;^" 1 / 2 ^)) 3 ) and g t G L x (0, +oo; (F" 1 / 2 ^)) 3 ) , (2.H11) 
/i G L°°(0, +oo; V') n L 1 (0, +oo; il) and h t G L 1 (0, +oo; V') . (2.H12) 

The above assumptions yield in particular 

F G L°°(0,+oo;W') and P t G 1^(0, +oo; W) , (2.38) 
h G L 2 (0, +co; V') . (2.39) 

Furthermore, using (|2.38p . it is not difficult to prove (see Remark 2.3] for all details) that 

3 Fqo G W : F(t) -» Foo in W as t -» +oo. (2.40) 

Hence, Theorem 12.11 ensures that for every quadruple of initial data ($o> u o> ^o) fulfill- 
ing (|2.1ip - (|2.14p there exists (at least) a solution trajectory i? s , u, X) : (0, +oo) — > F x 
-ff 1 (r c ) x W x i/ 1 (r c ) originating from (i?o> i?s) u o 5 ^o)- The ensuing Proposition 12.51 contains 
some suitable large-time a priori estimates for such trajectories. Such bounds shall enable us 
to conclude that the associated w-limit set (|2.44p is non-empty, and that its elements solve the 
stationary system associated with Problem (P) (see Theorem 12.21) . 

As we shall see, these results in fact hold for a class of solutions of Problem (P), namely 
approximable solutions which, in order not to overburden the paper, we shall precisely define in 
Section U] only (cf. Definition 14. lip . Here, we may just mention that the notion of approximable 
solution is tightly linked to the approximation procedure developed in Sec.[Uto prove the global 
existence of solutions to Problem (P) (see Theorem I2.ip . Such a solution notion allows us 
to perform rigourously on system (I2.27p - (|2.34p some of the a priori estimates on which our 
large-time analysis relies (cf. Remark 12.61) . 

Long-time a priori estimates. 

Proposition 2.5. Assume (|2JIT]) - (|2JI5|) and fTH9l) - (l2.H12p . Let (i? , #2, u , X ) be a quadru- 
ple of initial data complying with (|2.1ip - (|2.14|) . Then, there exists a constant K\ > 0, only 
depending on the functions A, k, a, and on the quantity 

M := IKollii(n) + lKllii(r c ) + lluollw + S(u ) + ||XolUi ( r c ) + P(^o)||ii(r c ) 

+ I|F||l°°(o,+oo;W) + ||Ft||Ll(o,+oo;W') (2-41) 
+ INlL°°(0,+oo;V)riL 1 (0 1 +° o ;- H ') + \\ht ||il(0,+oo;V) i 

such that for every approximable solution (in the sense of Definition ^. 17\ ) (i?, i? s , u, X) to Prob- 
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lem (P), fulfilling initial conditions (|2.23p - (|2.26p . there holds 

l|Vtf|| L 2(o,+oc;//) + l|Vtf s || L 2(o,+oo;L2(r c )) < #1 , (2.42a) 

II^IU^o.+ooiiHfi)) + \Ws\\l^ {o,+oo;^(r c )) < ^i ! (2.42b) 

||0-0-lb(o,+oo ; L>(r c )) <^i> (2.42c) 

ll u i||L 2 (0,+oo;W) + II u IIl°°(0,+oo;W) + ll«( u ) Hl°°(0,+oo) + II *i II L 2 (0,+oo;L 2 (r c )) < ^1 , (2.42d) 

ll x IU oo (o,+oo;fl' 1 (r c )) + ll/3(^)llL=°(o,+oo;L 1 (r c )) < -^l ) (2.42e) 

l|£Un(tf)|U2 (0j+oo; y + ||a t ln(^)|| i2(0)+00 ^i(r c )o < #1 ■ (2.42f) 



Furthermore, for all 5 > i/iere exist constants K2(S), K^(5,p) > 0, depending on 5, on the 
functions X, k, a, and on the quantity M (I2,4ip (K^(5, p) on p £ (0,2) as well), but independent 
of ?9 s ,u, X), such that the following estimates hold 

\W\\l™{8,+oc;V) + ll^s||L° c (5,+oo;// 1 (r c )) + IMlL°°(<5,+oo;ff 2 (r c )) , n ,„ . 

(2.43a) 

+ ||^t||L 2 (5,+t»;fl' 1 (r c ))nI ( ao (i5,+oo;L a (r c )) + ll u ll W^fr+oofN) ^ ^2^), 

\\Mms,+oo;L^/y { Q)) + ||^ s ||i2(5,+oo;L2-p(r c )) < #3(<M • (2.43b) 

Remark 2.6. In Section f3. II we shall give a formal proof of the above result, in which all the esti- 
mates leading to (j!T4"2"j) and (fiP3"j) shall be performed on the PDE system ([2^27]) - ([234]) directly. 
In particular, this shall involve the formal differentiation of equations (|2.29p and (|2.31|) . as well 
as formally testing (|2.27p . (|2.28p by $t an d c\i? s , respectively. All of these calculations shall be 
rigorously justified, by working on a suitable approximation of Problem (P), in Section [4.31 



Results on the w-limit of solution trajectories. Now, for a given quadruple of initial 
data (#o, $s, u o> Xo), complying with (|2.11|) - (|2.14|) . let ($, $ s , u, X) be an approximable solution 
starting from ($o, Uo, Xo), (its existence is ensured by Theorem 12. ip . We aim to investigate 
the cluster points for large times of the trajectory of # s , u, X), in the topology of the space 
H 1 " 6 ^) x iJ 1-e (r c ) x (H 1 " 6 ^)) 3 x i? 2 ~ £ (r c ), with an arbitrary e > 0. To this aim, we define 
the w-limit set u;(i?, # s , u, X) of the trajectory ($(i), $ s (i), u(i), X(i)) t >o a s 

w(0,0 8 ,u,X) :={(tfoo ) tf,,« ) Uoo I X 00 ) e y x i/ 1 (r c ) x W x ff^rj : 
3{t n } C [0, +oo), t n y +oo as n I oo, 

with (tf(i n ), i9 s (i„), u(i n ), X(t„)) -> (#oo, # s ,oo, u^, Xoo) *' 2 ' 44 ' ) 

in f 1 -^) x # 1_e (r c ) x (# 1 ~ e ($7)) 3 x # 2 - f (r c )}. 

For simplicity, we choose to omit the dependence on the initial data (i?o, u Cb Xo) (and on 
the parameter e), in the notation cj($, $ s , u, X). Notice that the latter would be replaced by 
the more customary lj^q, uo, Xo) if we additionally disposed of a uniqueness result for the 
approximable solutions to Problem (P). 

The following theorem shall be proved in Section 13.21 

Theorem 2.2. Assume ^MH\ - (l2~H5l) and (l2~H9l) - ()2.H12p . 

Then, for every quadruple of initial data ($o> uo, Xo) complying with (I2.11|) - (|2,14p and 
/or ewery approximable solution ($(t), , & s (t),u(t),X(t))t>o originating from (■&q,t)^,Uq,Xq), for 
all e > i/te associated uj -limit set i? s , u, X) is a non-empty, compact and connected subset 

ofH 1 -^) x ff 1 - 6 ^) x (ff 1 - 6 ^)) 3 x F 2 ~ e (r c ). 
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Furthermore, every quadruple ($oo, $s,ooi u oo> ^oo) £ i? s , u, X) is a solution of the sta- 
tionary system associated with Problem (P), namely 



(2.45b) 



Xoo G # 2 (T C ) and < 



3 $00 > : $00(2:) = "#oo /or o-e. x € 0, i? s ,oo(^) = $00 /or a.e. x G r c; (2.45a) 

a( Uoo ,v) +tfoo / div(v)+ / (X+ Uoo + 77^) • v = Foo VveW, 
Jn Jr c 

Voo G a( Uoo ) zn (tf" 1 / 2 ^)) 3 , 

-AXoo + £oo + o-'(Xoo) = -A'(X 00 )'!9 00 - ^Cooluool 2 a-e. in T c , 
£00 £ ^ 2 (r c ), Coo £ /?(X 
Coo £ ^°°(r c ), Coo £ tf(Xoo) a-e. in T c , 
c?n a ^oo = a.e. in dT c . 

(2.45c) 

Remark 2.7. Let us emphasize that, in the limit as t — > 00, the system is in a state of 
thermomechanical equilibrium. Indeed, the dissipation, described by the pseudo-potentials (II. 3D 
and (fL~4"j) . has vanished in (|2.45ap ~ (|2.45c|) . 

No uniqueness result is available on the stationary system (|2.45a|) - (|2.45c|) . Hence, one cannot 
deduce directly from Theorem 12.21 that u;(i?, $ s , u, X) is a singleton and that the whole solution 
trajectory (i?, # s ,u, X) thus converges, as £ — > +00, to a unique equilibrium. However, the next 
result (whose proof is postponed to Section l3.2p shows that, under more specific assumptions 
on the operator (3 and on the nonlinearities A and cx, it is possible to uniquely determine the 
X-component of the elements in # s , u, X). 

Corollary 2.1. Under assumptions (j2.Hlh n2.H5j) and (j2.H10j) - (j2.H12j) . suppose further that 

(3 = 91r m m »i, for some —00 < to* < to* < +00, (2.46a) 
A is non- decreasing on [to*,to*] ; (2.46b) 
a'(x) > for all x G [m*,m*]. (2.46c) 

Then, for any quadruple (i?oo, $s,oo, Uoo, ^00) £ w(t?, $ s , u, X) t/iere /10/ds 

X 00 (x) = m, Vier c , (2.47) 

and we have as t — > +00 

X(i)->Xoo in H 2 - e (T c ) for all e> 0. 

Remark 2.8. Corollary 12.11 ensures that, in the case when the latent heat is positive, (3 has 
a bounded domain (which is the interesting case from a physical point of view), and cohesion 
in the material (which is included in the decreasing part of a, see (jl.2p ) is not too large with 
respect to the remaining part of the potential a, the glue tends to be completely damaged in 
the large-time limit. We remark that this is the result one would expect from experience. 

3 Proofs 

Notation 3.1. Henceforth, for the sake of notational simplicity, we shall write (-, •) for all the 
duality pairings W /(-, -) w , y,{-, -) v , and #-i(r c )'("j ')h 1 (f c )^ anc ^ further, denote by the symbols 

c, c' C, C most of the (positive) constants occurring in calculations and estimates. 

Notice that the above constants shall not depend on the final time T. 
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3.1 Proof of Proposition 12.51 

The following is a formal proof, based on two main a priori estimates on system (|2.27p - (|2,34l) , 
which shall be revisited in Section 14.31 Therein, we shall also rigorously prove the further 
solution regularity ()2.36ap n2.36ep . 



First (formal) estimate. We test (12371) by (l2~28l) by tf s , (l2~29l) by u t and (|2~3Tj) by X t , 

add the resulting relations and integrate them on the interval (0, t), with t € (0, +oo). Now, we 
take into account the formal identities 

*(^lnW^) = ||^)||Li(n)-||^o||L 1( n), 

(3.1) 

{dtH0'),#s) = \\MmiA<r e )-W\v<r e ), 

and the chain rule for the convex functionals 2, and (5 (cf. with |!Uj Lemma 4.1] and 
Lemma 3.3], respectively), and for the smooth function <r, yielding 



f 77 • u t = S(u(t)) - S(u ) , [ [ £x t = [ (3{X{t)) - I P{Xq) 
o Jr c Jo Jr c Jt c Jt c 

f a'(X)X t = [ a(X{t))- [ a(X ). 
io Jv c Jv c Jr c 

In the same way, an integration by parts and the chain rule for p(-) = (-) + lead to 
X+uu t = - f [ X + d t \u\ 2 



(3.2) 



o Jr c z Jo Jr. 



(3.3) 



where we recall that £ € H(X) = <9p(X). Finally, we observe that, by the properties ([2 
and (|2.5p of the forms a and 6, respectively, we have 

ra(u,ut)>^||u(t)||^-ilir a ||uo||^, Akih) >C b Tllutllw- (3-4) 

JO 1 1 JO JO 

Collecting (I3,ip ~ (l3,4p and observing that some terms cancel out, we get 

l|0(t)IUi ( n) + 1 W^Wh + J o J r k(x){# - tf s ) 2 + \\Mt)\\mr c) + J q l|w s ||| 2(rc) 

+ C 6 jT \\n t f w + ^\\u(t)\\ 2 w + 3(u(*)) + \J X + {t)\u{t)\ 2 



+ j Q ll^lli2 (rc) + 2l|VX(t)||i 2(rc) + ^ p(X(t))+(7(X(t))) 
< ll^olki(n) + ||tf°lki(r c ) + ^Kllw + «(uo) + ^||VX ||£ 2(rc) 



+c||Xo||i2 (rc) ||uo||w + ll^(^o)[Ui(r c ) + h(X )\\^(T c ) + h + h , 
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where have also used the continuous embedding (|2.2p for the third term on the right-hand side 
of (|3.3|) . Then, we estimate 

<M>< / \\h\\ v ,\\# - m(#)\\ v + f \\h\\v\\m(0)\\v 
o Jo Jo /o «\ 

rt i rt 1 rt ( 3 - 6 ) 



<cJ o \\hfv> + ll l|W|l ^ + |^i WHvWIliW), 



the second inequality due to the Poincare inequality for functions with zero mean value, and, 
with an integration by parts, we get 



h = ( (F,u t } = - / (F t ,u) + (F(t),u(i)}-(F(0),u ) 
Jo Jo 

/t (J I 

||F t ||w'|l u llw + -^||u(*)|lw + 2 H U ° Hw + C|l F lli«>(0,+oo;W') • 



(3.7) 



Hence, taking into account our assumptions on the data (|2.H10l) - (12.H12p (which yield (12351) 
and ()2.39p ). as well as (I2.1ip - (|2.14p . we may apply a variant of the Gronwall Lemma (cf. [HI 
Lemma A. 5]) and conclude estimates (|'2.42a|) . (|2.42bp . and ()2.42dp . while (|2.42cp follows from 
the bound for (fe(X)) 1 / 2 (t? — $ s ) in L 2 (0, +oo; L 2 (T C )) and the fact that k is bounded from below 
by a strictly positive constant, see (|2.H5p , Finally, thanks to (|2.H9p and (|2.9ap . we deduce 
from (ESI) that 



3C,c,c'>0 Vt€(0,+oo) : c\\X(t)\\ 2 L2{Vc) - c' < J (j3(X(t)) + ff(X(t))) < C . (3.8) 

Joint with the bound for VX in L°°(0, +oo; L 2 (F C )), this yields the estimate for ||X||^oo(o i + o;_ff 1 (r c ))- 
A fortiori, in view of (I2.H9P we also recover the bound for f3(X), and H2.42e[) ensues. In the end, 
estimate (|2.42f|) for e?tln(i?) and for 9fln(i? s ) follows from a comparison in equations (|2.27l) 
and (I2.28|) . respectively. For example, using (|2.H5|) and (|2.42cj) . (I2.42el) one observes that 

\\k(X)($-$ s )\\ L 2 {0t+oo . vl) < C(||X|| L ^ (0i+oo;L 4 (rc)) + 1)||^ -^|| L2(0i+oo;L2(rc)) < C . 

Hence, in view of the bound (I2.42al) for W in L 2 (0, +cx>; H), of (|2.42dl) for div(u<) in L 2 (0, +oo; H 3 ), 
and of (|2.39p . one concludes the estimate for dt ln(t?) in L 2 (0, +co; V). 



Second (formal) estimate. In what follows, we shall formally treat the maximal monotone 
operators a, H, and (5 as nondecreasing and Lipschitz functions. Indeed, the following estimates 
can be rigorously justified as in Section O regularizing these nonlinearities by their Yosida ap- 
proximations (which are in fact Lipschitz functions). Further, we shall work with a solution 
(f?,i? S ) u, X) enjoying the further regularity (|2.36p (see also (I2.37P ). Formally proceeding, for all 
5 > we let a constant K$, only depending on A, k, a, on the quantity M (|2.4ip . and possibly 
on 5, such that 

\\m\\v + \Ws(S)\\m { r c ) + KWIIw + \\Uml H r c) < Ks ■ (3.9) 

This formal assumption shall be discarded once we put forth the rigorous arguments of Section 
HI see Remark I3.2I for further comments. Hence, we are in the position of performing the 
following calculations. We test (l2~2TD by °d t and RTM by d t "9 s , differentiate RTM w.r.t. time 
and test it by uj, and differentiate (|2.3ip w.r.t. time and multiply it by X%. We add the 
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resulting relations and integrate them on the interval (5,t), with t € (<5, +00), also adding 
1/2(11^^)11^!^^ + ||f?s(t)||^i( r \) to both sides. Indeed, the Poincare inequality for the zero mean 
value functions yields 



c P \\m\\l <\(\\vm\\ 2 H + \\m\\h m 

C P \\Ut)\\ 2 H H T c ) < \ (\\VMt)\\h(r c ) + \\Ut)\\l H r c ) 



(3.10) 



for some positive constant Cp independent of t 6 (0, +00). We also notice that, for some other 
constant c also depending on the embeddings (|2.2f) (|2.3H . there holds 

c \m - Mt)\\l H T c) < % (iw*)u* + n^)ii Wc)) for a11 1 G ( ' +o °)- ( 3 - n ) 

Further, we remark that 

f [ d t (X+u)u t = [ [ X+\u t \ 2 + [ [ H(X)X t uu t . (3.12) 
Js Jr c Js Jr c Js Jt c 

Taking into account the cancellation of some terms and the coercivity and continuity of the 
forms a and b (|2.4|) - (|2.5|) . using the formal identities 

ft ft f |„n |2 ft 



L 



(dM^t) = J s JJ-^f, J s (dtM#s)A*8) = J s j v (3.13) 



as well as (|3.10p - (|3.12p . we end up with 



(3.14) 



s Jn # 2 11 J s J Tc v /v °' tv J s Jr c #s 

+ ^¥s(t)f m{Tc) + c\\m - Ut)\\l H r c) + y K(*)Hw + C a f 

+ [ [ X + \u t \ 2 + [ (a'(u)u t ,u t } + hx t (t)\\ 2 L2{rc) + [ [ \VX t \ 2 

+ f [ P'(X)\X t \ 2 < CK 5 + K 2 + h + h + h + h + h + h , 
Js Jr c 

in which we have controlled the term l/2(||i?(i)||| 1 ^ + ||^ s (t)||^i( r )) on the right-hand side 
by (|2.42bj) . Integrating by parts, we have 



<Mt> = ~ J (ht,#) + (h(t),m) - (h(S),HS)) 

f ||MHh% + CINll^o.+oo^) + ^¥{t)\\v + K S , 



(3.15) 



while we estimate 



I 



1 



f (F t ,u t ) < [ ||F t || W '||ut||w, (3.16) 
Js J 8 

-2 H(X)X t uu i < 2 / ||X t || i 2( rc )||u|| £ ,4( rc )||u t || L 4( rc ) 

Js Jr c Js 

< C|I u IIl°°(o,+oo ; w) yj^ \\Xt\\h t (r c )+ J g J r ll u *llw^ <c, 

-\f s j v H'(X)|X t | 2 |u| 2 <0, (3.18) 



(3.17) 
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(|3.17p due to the continuous embedding (|2.2p and to estimate (|2.42dj) . and (|3. 18j) following 
by monotonicity. Moreover, recalling (|2.H3p . (|2.H4p . and ()2.42dp . and applying the Holder 
inequality 



h = - f [ a"(X)\X t \ 2 <L a Kf 
J8 JT C 



\" (X)\X t \ 2r d s < L\ [ IIXtllia^jllXtll^crcjIl^llL^c) 
r c Js 

2 



< - 

~ 4 



lj s ll vx tlli 2 (r c ) + -f + C j 5 H Xt Hi 2 (r c )ll^llHi(r c ) 



(3.19) 
(3.20) 



C also depending on the constant of the continuous embedding (|2.3|) . Finally, we integrate by 
parts the third term on the left-hand side of (13.141) . so that 

k(X)(0 -•&.)%(•&-•&.) 

r c 

= \J t k{x(t))\m-ut)\ 2 -\f r k(x(s))\m - us)\ 2 - h 

> f IW*) - Mt)\\l H r c) ~ c (\\X(8)\\ LHrc) + l) (\\m\\ 2 v + ll^)ll^(r c) ) ~ h 

> f ¥{t) ~ Ut)\\% { T c ) -C(K S + 1)-I 9 , 

the latter inequality due (l2TH5j) . to (l23cj) . (|2.42e|> . and (USD, while, using (|2.H5|) and (|2.42dll . 

we estimate 



= \ \ j 



1 

2 

< — 
~ 2 



\MlHt c )¥ - MlHf c )\W - MLHr c ) ( 3 - 22 ) 



4/l|vx t ||| 2(rc) + f + c£ 



^sllL2(r c )ll^ ~~ ^sllz,4(r c ) 



Next, we collect (I3.14p - (l3.22p . observing that the ninth term on the left-hand side of (13.141) 
is non-negative, and so are the tenth and the thirteenth terms, by monotonicity of a and f3, 
respectively. Then, taking into account the summability properties (I2.38H and (|2.39l) . as well 
as estimate ()2.42dp for \\X t ||L2(o,+oo ; L2(r c )) an d (j2.42cp for - ■& s \\l2(o,+oo;L^(t c )), we apply 
the standard Gronwall lemma and its abovementioned variant (cf. (9j Lemmas A. 3, A. 5]) and 
conclude estimate (|2.43ap for # s , u, and Xt. A comparison in (|2.3ip yields (possibly for a 
larger 1(2(8)) 

\\-AX + al^s, + oo;mr c ))<K2(8), 

whence an estimate both for ||CllL oo (5,+oo;L 2 (r c )) an d for || — AX\\ Lao rg t+Q0 . L 2rr c \\ °y the monotonic- 
ity of f3. Then, the bound for ||X|| Loo (s,+oo;H 2 (r c )) follows from elliptic regularity. Furthermore, 
the estimate on the first and on the fourth integral terms on the left-hand side of (|3.14p gives 
(notice that $ > a.e. in Q x (0, +00) and $ s > a.e. in T c x (0, +00)) 

\\dt$ 1/2 \\m8, + oc;H) + l|9^ s 1/2 ||L 2 ( 5 , + oo;L 2 (r c) ) < C(8) . (3.23) 

On the other hand, due to the continuous embeddings V C L P (Q) for all 1 < p < 6 and 
^(Tc) C L q (T c ) for all 1 < q < oo, (|2.43ap in particular yields 

VI < q < oo 3C q (8) > : \W 1/2 \\ L °°(s,+oo;L^(n)) + ll^ /2 ||L-(5,+oo;L 9 (r c) ) < C q (8) . (3.24) 
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Hence, (|2.43bj) follows from (|2.43a|) . (f3T23j) . (f334"|) . and the Holder inequality, giving 

WMl^q) < W II2 \\h¥ 1/2 \\l^) 

WWbWv-pFc) ^ H^^ /2 |lL 2 (r c )ll^ /2 |lL(4- 2 p)/P ( r c ) f o r all < p < 2. 
This concludes the proof. 

Remark 3.2. Following [B], we point out that a possible way to perform estimate (|3.14|) more 
rigorously, without assuming (|3.9p . would be to fix a smooth "cut-off" function q : [0, +oo) — ► 
[0,+oo), with q' G L°°(0,+oo) (for example, q(t) = tanh(t) for all t > 0), and test (l2~271) by 
(|2.28|) by qdt"& s , the time derivative of (I2.29P by <nit, and the time derivative of ()2.3ip by 
qX-t. However, to keep calculations simpler we have postponed this procedure to the rigorous 
proof of Proposition 12.51 in Sec. HI 



(3.25) 
□ 



3.2 Proof of Theorem [Ml 

In view of estimates (12.42dp and (I2.43ap . for all 6 > the trajectory {(#(t),# s (t),u(i),X(i)),t > 
5} of every approximable solution ($, $ s , u, X), starting from a quadruple of initial data ($o> ^! uo, Xq) 
as in (|2.1ip - (|2.14p . is bounded in V x i/ 1 (r c ) x W x H 2 (T C ). Hence, it is relatively compact in 
H 1 " 6 ^) x H l - e (T c ) x (H l - e (tt)Y x iJ 2_e (r c ) for all e > 0, and with a standard argument one 
finds that the set $ s , u, X) is a non-empty and compact subset of the latter product space. 
Thanks to flZgfr , w(0,0 a ,u,X) is connected in H 1 ^^) x fT 1-6 ^) x (tf 1 "^)) 3 x F 2 - £ (r c ) 
for all e > as well, by a well-known result in the theory of dynamical systems, see e.g. j!5j . 

Now, let us fix (i?oo, $ s ,oo> Uqo, Xqo) € u, X) and an increasing sequence {t n } C 

(0, +oo) such that t n /" +oo as n — > oo and 

<# s (7j n ),ll(f n ),X(7j n )) -> (^oo^s.oo.Uoo^oo) 

in H 1 -'^) x (r c ) x (if 1 "^)) 3 x H 2 - € (T C ). 

Following a well-established procedure, for all n G N and T > we introduce the translated 
functions 

d n {t) :=d{t + t n ), $ Stn {t):=d s {t + t n ), Un(t) :=u(t + t n ), X n (t) := X(t + t n ) , 
for t G [0,T]. We also set for a.e. t G (0,T) 

r7 n (t) := rj(t + t n ), Cn(t) := £(i + *n), £»(*) := C(* + *n), F n (i) := F(t + t n ), /t n (rj) := fr(i + rj n ). 

Clearly, r? n (rj) G a(u n (i)) for a.e. i G (0,T), and £„(x,i) G (3(X n (x,t)), ( n (x,t) G H(X n ( s,t)) fo r 
a.e. (x,t) G T c x (0, T). For later convenience, we point out that, in view of (|2.38|) and ()2.H12p . 
for all n G N and T > 

l|Fn||z,°°(0,T;W") < ||F II L°° (0,+oo;W) > (3.27) 
||^n|U°°(0,T;V) < |NU°°(0,+oo;V)> 1 1 ^ II (0,T;/f) < II ^ IIl 1 (0,+oo;fT) ■ (3.28) 

Furthermore, (12.391) leads to 

T rt n +T r+oo 

\\h n \\ 2 v ,= / \\h\\ 2 v , < / \\h\\ 2 v , ^0 as 7i — ► co (3.29) 
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whereby 

h n -> in L 2 {0,T;V), h n ^* in L°°(0, T; V'), (3.30) 
whereas, using (|2.40[) . we verify that 

3 Foo € W : F n (t) -» Foo in L p (0, T; W) for all 1 < p < oo . (3.31) 

Now, 

s,ni Uji, ^ri) ^7ru £ni Cn) is a solution to Problem (P) (with data h n and F n in place of 
h and F) on the interval (0, T), and, by construction, the functions # n , ?9 Sjn , u n and X n comply 
the initial conditions 

tf„(0)=tf(t n ), u n (0) = u(t n ) inO, 0., n $)=Mtn), X n (0)=X(t n ) mT c . (3.32) 

Exploiting (|3.26p . we shall prove that the quadruple (i?ooi $s,oo, u ooj ^oo) fulfils the stationary 
system (|2.45p by passing to the limit as n — ► oo in the abovementioned initial-boundary problem. 

To this aim, in view of estimates (|2.42p and (|2.43p (indeed, we may suppose without loss 
of generality that, e.g., t n > 1 and choose p = -y in (|2.43bp ). we remark that there exists a 
constant > 0, independent of n E N and T > 0, such that 

ll$ra||L°°(0,T;V) + || V# n ||l2( 0) T;H) + ll^t ln(l? n ) \\l^(0,T;V) + l|3fc$n llL2( 0i T;L 12 /7(n)) ^ K 4 , 
\Ws,n\\L°°(p,T;H 1 (J' c )) + ll V ^,n||L 2 (0,T;L 2 (r c )) < K 4 , 
\\9t'^-('&s,n)\\L'(p,T;H 1 (r c y) + ll^^s,n|li 2 (0,T;i 12 / 7 (r c )) ^ ^4, 

||#n ~ #s,n||L 2 (o,T;L 2 (r c )) < K4 , (3.33) 

||u„||vi/l>°°(0,T;W) + ll^t u n||L 2 (0,T;W) + II «( u n) ||l°° (0,T) < ^4 , 

ll^n||L°°(0,T; J H' 2 (r c ))nVl/ 1 '° o (0,T;L 2 (r c ) + ll^t^n Hi 2 (CT;^ 1 (r c )) + ll/?(^n) ||L°°(0,T;Z, 1 (r c )) 

+ ll£n||L°°(0,T;L 2 (r c )) ^ ^4 • 

Moreover, since < £ n < 1 a.e. in T c x (0, T) there holds 

llCnlU°°(r c x(o,T)) < 1 ■ (3.34) 

Finally, let us prove a further estimate independent of n € N but depending on T > 0. Exploit- 
ing (I3.33|) (which also yields a bound, independent of n S N and T > 0, for ||X+u n || LOO ( r.#-i/2(r c ))) 
and (|3.27|) . we argue by comparison in (|2.29|) and conclude that 

\\Vn\\mo,T;H-V2(r c )) < K 5 (T) for all n £ N. (3.35) 

Estimates (I3.33p - (I3.35H . joint with standard weak-compactness arguments, the compactness 
results \18\ Thm. 4, Cor. 5], and the Ascoli-Arzela theorem in the framework of the weak 
topologies of W and ^T 1 (r c ), yield that there exist a subsequence of {(# n , $s,n, u n , X n , rj n , £ n , ( n )} 
(which we do not relabel), and functions ?? s , u, X, f), £, Q for which the following convergences 
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hold as n — > oo 

tf n -*tf mL^(0,T;V)nH 1 (0,T;L 12 / 7 (n)), 

n i (3.36) 

$n — * $ in &([0,T\;H 1 - p (Sl)) for all p G (0,1), 

tf s ,„^*tf s inL oo (0,r; J ff 1 (r c ))n J ff 1 (0,T;L 12 / 7 (r c )), 

tf s , n ^ in C°{%T]-H l ~P{T c )) for all p G (0,1), 

u n ^u in H l (Q,T; W) 

u n ^u in C°([0,r]; J H" 1 -''(fi) 3 ) for all p G (0, 1) , 

rj n ^fj inL^O,^^- 1 / 2 ^)) 3 ), (3.39) 

x n ^*x inL 00 (o,r;F 2 (r c ))nF 1 (o,r; J ff 1 (r c ))n^ 1 ' 00 (o,T;L 2 (r c )), 

X n ^X in C°([0,T};H 2 -P(T c )) for all p G (0, 2), 

£ n -*£ inL°°(0,r;L 2 (r c )), (3.41) 
Cn^* C in L°°(0, T; L p (r c )) for all 1 < p < oo . (3.42) 

Clearly, from the positivity of the sequences {i? n } and {$ s ,n} we deduce that 

0(x,t)>O for a.e. (x,t) G O x (0,T), i9 s (x,t)>0 for a.e. (a?, t) G T c X (0, T) . 

Furthermore, arguing in the same way as for (|3.29|) - (|3.30|) . one sees that estimate (|2.42f|) for 

||9tln(i?)||£,2( 0j+oo .y/) and \\d t ln(^ S) „) Hl^o.+oo;^ 1 ^)') entails that, as n -> oo, 

d t ln(tf n )^0 inL 2 (0,T;y'), d t hi(tf s , n )^0 in L 2 (0, T; i? 1 (r c ) / ). (3.43) 
Analogously, (|2.42dp implies that 

d t u n ^0 inL 2 (0,T;W), (3.44) 

and 

S t X n ^0 inL 2 (0,T;L 2 (r c )). (3.45) 
In the same way, in view of (12.43b|) for $t and dt"& s , we conclude that 

d t $ n ^0 inL 2 (0,T;L 12 / 7 (f))), dtf., n -» in L 2 (0, T; L 12 / 7 (r c )). (3.46) 

Finally, from (|2.42ap and (|2.42c|) . we infer that 

(0»-0»,„)-O inL 2 (0,T;L 2 (r c )) (3.47) 
W„-»0 mL 2 (0,T;H), W«,„ -»> in L 2 (0, T; L 2 (r c )). 

With (|3.44|) we get lit = a.e. in 0, so that u is constant in time in f2. Hence, 

u(x,t) = u(x,0) = lim u(x,i n ) = u 00 (x) for a.e. x G for all t G [0, T] 

n— >oo 

and analogously, from (|3.45|) and (|3.46p we deduce for all t G [0,T], 

$(a;, i) = $oo(x) for a.e. rc G 0, , $ s (x, t) = "& s ,oc(%) ) = ^oo(a^) for a.e. x G T c . 

Finally, owing to (|3.4T|> we get that #oo and i? Sj00 are constant (as V#oo = V# Sj00 = 0) and 

^oo| r = ^s,oo a.e. inr c . (3.48) 
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We also point out that (|3.38p and (|3.42p yield, as n — > oo, 

Cn\u n \ 2 ^* C|Q| 2 in L°°(0, T; L 2 ~ P (T C )) for all < p < 2. (3.49) 

Further, combining (|3.46j) with (|3.33p and (|2.H4|) . we easily conclude that dtX(X n ) = A'(X n )i9 t X n — 
in L 2 (0, T; L 2 ~ P (T C )) for all p £ (0,2). Arguing in the same way as in the proof of [U 
Prop. 4.4], we use (I336l) - (l3^0l and (I3^3l) - (l09l to pass to the limit in Problem (P) and, 
taking into account (|3.30p - (|3.3ip . we conclude that, almost everywhere in (0, T), the functions 
($oo,#s,oo,Uoo, Xoc^^C) fulfil fl2.45a.H and 



a( Uoo ,v)+ / t? 00 div(v)+ / (X+ Uoo + ^) • v = Foo VvGW, 



(3.50) 



- AXoo + £ + ^(Xoo) = -A'^oo^oo - -C|uoo| 2 a.e. in T c , (3.51) 

joint with the no- flux boundary conditions (|2.34l) . Combining convergences (I3.40p . ()3.4ip . 
and (|3.42p with the strong-weak closedness of the graphs (3 and H, we find that 

Ce^Xoo), CeH(X 00 ) a.e. inr c x (0,T), (3.52) 

while, in view of the maximal monotonicity of the operator (induced by) a on L 2 (0, T; (H~ 1 / 2 (T C )) 3 ) 
and of [H Lemma 1.3, p. 42], we deduce that 

f7Ga( Uoo ) in {H~ l / 2 {T c )f a.e. in(0,T) (3.53) 

from the inequality 

limsup / / rj n u n < / f\ ■ vl^ for alii G (0, T) , 

?woo Jo Vr c ./o Vr c 

which can be verified in the same way as in the proof of [U Prop. 4.4]. 

Remark 3.3. Notice that (i?oo> "$s,oo 5 Uoo, X^) solve the stationary problem ()2.45a|] - ()2.45cp as- 
sociated with the evolution system (|2.27p - (|2.34p . Moreover, we observe that (|3.44fl — (|3.47j) entail 
that dissipation vanishes in the limit as t — > oo (see f 1 1 . 3 [) and (II. 4p ). 



Proof of Corollary 12.11 For the sake of completeness, here we repeat the same argument 
developed in the proof of [31 Prop. 2.5]. We test the first of (|2.45cp by (X^ — to*) and integrate 
on T r . We obtain 



(Xoo - m*) < , 



|V(Xoo-m*)| + / £00(^00 - m* 
'r c 



jf ^'(Xoo) + A'(X 00 )^ i00 + iCoo|u o| 2 



the latter inequality due to rt2.46b|) — (j2.46cH and the fact that to,* < X^ < to* a.e. in r c . On 
the other hand, the second term on the left-hand side of the above inequality is non-negative 
by monotonicity, so that we deduce that V(Xoo — m,) = a.e. in T c . Thus, there exists 
some constant g > m, such that = g a.e. in T c . Now, integrating (|2.45cj) and again 
recalling (I2.46bp - (l2~l6cll . we find that 



Coo < 0. 

Hence, necessarily g = ?n*, and ()2.47p ensues. □ 
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4 Rigorous estimates 

This section is devoted to the rigorous proof of Theorem 12.11 and of Proposition 12, 51 Thus, 
in Sec. 14.21 we shall specify the variational problem (depending on two parameters e > and 
fi > 0) approximating Problem (P) and outline the main steps of the proof of its global well- 
posedness. Hence, in Sec. 14. 31 we shall prove some estimates on the approximate solutions, which 
are independent of the parameters s and fx and in fact hold on (0, +oo). This shall enable us 
to pass to the limit in the approximate problem first as e \ 0, and secondly as /i \ 0, and to 
conclude the proof of Theorem 12.11 in Sec. 14.41 and of Proposition 12.51 in Sec. 14.51 

Since the (double) approximation procedure for Problem (P) strongly relies on the usage of 
Yosida regularizations of the nonlinear operators In, a , /3, and of the Heaviside operator H, in 
the following section we recapitulate some related preparatory results. 



4.1 Recaps on Yosida regularizations 

Regularization of In. For fixed fi > 0, we denote by 

: = (Id + fxln)" 1 : R^R (4.1) 

the resolvent operator associated with the logarithm In (where Id : R — > R is the identity 
function), and recall that r M : R — > (0, +oo) is a contraction. With easy calculations, one 
sees that r^(l) = 1, so that, by contractivity, there holds 

i>(x) < |x| + 2 for all x € R. (4.2) 

The Yosida regularization of In is then defined by 

Td — r 

ln M := — £ : R -> R . (4.3) 

It follows from [9l Prop. 2.6] that for all \x > the function ln At : R — > R is non-decreasing 
and Lipschitz continuous, with Lipschitz constant l//i. 

For later convenience, as in [5] we also introduce the following function 

^(x):= f X sln'^ds. (4.4) 
J o 

We point out that, since is decreasing on (— oo,0) and increasing on (0,+oo), there 
holds for all fi > 

^(x) > ^(0) = for all x G E. (4.5) 

The following result collects some properties of ln^ and ^ which shall play a crucial role 
in the proof of Proposition 14.51 

Lemma 4.1. The following inequalities hold: 

2 

3u*>0: V>e(0,u*) Vx>0 In' (x) < - , (4.6a) 

p x 

M ii > Vx G R In' (x) > — j— . (4.6b) 

MV ; ~ |x| + 2 + n y ' 

As a consequence, J?n satisfies 

3^>0: V/i£(0,/i«) Vx>0 ^ M (x)<2x, (4.7a) 
3d, C 2 >0: V/i>0VxeM JJx) > d\x\ - C 2 . (4.7b) 
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Proof. Ad (|4.6p . Using the definitions (|4.ip and f|4.3j) of r M and ln^ and repeating the 
calculations in the proof of [SJ Lemma 4.2], it is possible to show that 

In' (x) = — for all x G R, 

r M (rr)+ / u 

which, combined with (14. 2D . yields (|4.6bj) . For the proof of ()4.6ap . which follows the very 
same lines, we directly refer to Lemma 4.2]. 

Ad (|4.7|) . Estimate (|4.7a|) is an immediate consequence of (|4.6ap and of the definition of 
t/u- We shall now prove ()4.7bp for x > (the inequality in the case x < being completely 
analogous). Indeed, from the inequality 



sln'(s)>l fl+ foralls>0 



(which is an immediate consequence of (I4.6bp ). we deduce that 
S»{x) >x- f X M + 2 ds 



o s + ii + 2 

= x-(ji + 2) ln(x + (i + 2) + (/i + 2) ln(// + 2) > Cix - C 2 , 
for some suitable positive constants C\ and Ci- □ 

Regularization of a. For fixed \i > 0, we shall denote by 

a M : (iJ 1/2 (r c )) 3 -» (iJ _1 / 2 (r c )) 3 the /i-Yosida regularization of a, 

and recall that, by pQ Prop. II.l.l], the operator a M : (iT 1 / 2 ^)) 3 -» (H~^ 2 (T C )) 3 is 
single- valued, monotone, bounded and demi-continuous. Being a a non-negative functional 
by (!2.H2p . its primitive 

S„(u) := min_ ( — V ^ gl/2 ( r °) + S( v )^ for all u G (H 1 / 2 ^)) 3 

vefl(a) \ ^ J 

(which is Frechet differentiable on (i7 1 / 2 (r c )) 3 ), fulfils 



< S M (u) < S(u) for all u G D{a). (4.8) 

Regularization of /3. We shall also use 

(3u : R — > M the /i-Yosida regularization of /?. 

With straightforward calculations one verifies that, thanks to ()2.H9p . the Yosida approxi- 
mation 

3u(x) := min ( W-ZjEI + p(y)) for all x G R 



of /3 verifies 

(Vi2>0 3C R >0 V/i>0Vx>0 : /J„(x) > ^+T x2 ~ C R > (4 
| ]C>0 V/O0VK0: ^(x) > |£-C7. 
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Regularization of H. For fixed \x > 0, we shall denote by 

Hji = P« : R — > R the /it- Yosida regularization of H, 
: R — > R the ^-Yosida approximation of (-) + . 

In particular, it can be checked that 

'0 ifx<0, 

ii0<x<fi, (4.10) 



:r 2 



2u 

x — — if X > /X. 



Using the definition of fLj and p^, it is straightforward to verify that 

< H M (z) < 1 for all x G R, (4.11) 
< p M (x) < (x) + for all i£l, (4.12) 

4.2 Approximation of Problem (P) 

A double approximation procedure. Let e, fi > be two strictly positive parameters. We 
consider the approximation of Problem (P) obtained in the following way: 

1. we add to (|2.27p the regularizing viscosity term e3?(0j) and to (|2.28[) the viscosity term 
e3lr c (dt'& s ) (3£ and 3?r c being the Riesz operators introduced in Notation 12. ip ; 

2. we replace the operators a in (|2.29p . (3 and H in (|2.3ip with their Yosida regularization 
a M : (H l l 2 {T c )f -> (H-^ 2 (T C )) 3 , ^ : R -► R, and H M : R -» R; accordingly, we replace 
the term X + u in equation (I2.29P by p^(X)u; 

3. both in (|2.27)) and in (|2.28p we replace the logarithm In with its Yosida regularization ln^. 



Approximate initial data. In order to properly state our approximate problem, depending 
on the parameters e > and fi > 0, we shall need some enhanced regularity on the initial data 
for and S . The following result concerns the construction of sequences of suitable approximate 
initial data {0^} and {0g £ }, which in fact depend on the parameter e > only. 

Lemma 4.2. Assume that the initial data 0o and 0^ respectively comply with ([2. lip and (|2,12p . 

Then, 

1. there exists a sequence {0^} C V such that for all e > and fx € (0, /x*) ((/,* being as in 
Lemma \4-1\ ) there hold 

e 1/2 \W° £ \\v<\\M\v, (4.13a) 

SM)< 2 \\MlH{i), (4.13b) 

n 

i?° -» O *» ^(n) as e \ , (4.13c) 

ln M (0°) -» ln M (0 o ) in as e \ ; (4.13d) 
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2. there exists a sequence {p° s e \ C H l {T c ) such that for all e > and fi G (0, /i*) 

e 1/2 ||< £ ki(r c) <ll^ki(r c y, (4.14a) 
J^« £ )<2||<|| il(rc) , (4.14b) 



e #2 *n L 5 (r c ) as e \ , (4.14c) 
M<e) -f ln M (i?2) in L«(r c ) as e \ . (4.14d) 

Proof. We shall carry out the construction of the sequence {t?|?}, the passages for being 
completely analogous. For all e, \x > 0, we let $ £ € 1/ be the solution of the variational equation 

f + £l/2 I Vtf°Vv = / $ v for all v G F (4.15) 
Jn Jn Jn 

(notice that ft® € V is well-defined, thanks to (l2~T5l) ). Then, (|4.13aj) can be straightforwardly 
proved, together with 

-f i? in V" as e \ 0. (4.16) 
Furthermore, the maximum principle shows that, being i?q > a.e. in thanks to the second 

of (J2TU), 

> a.e. in Q. 

Then, testing (|4.15p by 1, one obtains 

II^IUi(O) = ll^olUi(Q) for all e > 0. (4.17) 

Therefore, (|4.13bp follows from combining (|4.17|) with (|4.7a|l . Finally, to check (|4.13cp we 
test (I4TT5D by (tf^" 1 and with easy calculations deduce that H^Hz^o) < ||^o||iP(n)- Com- 
bining this estimate with (|4.16|) using the uniform convexity of LP(fl), we conclude the strong 
convergence of {$®} in L^(fi). Hence, (I4.13dj) is a direct consequence of (I4.13dj) and the Lipschitz 
continuity of ln^. n 

Variational formulation of the approximate problem. We thus obtain the following 
boundary value problem, which we directly state on the half-line (0, +oo) in view of the long- 
time a priori estimates of Proposition 14.51 

Problem (Pe). Given a quadruple of initial data (i?^, E , Uo, Xo), uo and Xq being as in (|2.13l) 
fl2HU), and #° and i?° 6 fulfilling (j4"7L3l) and (|4TT4"j) respectively, find functions (i9, i? s , u, X), with 
the regularity 

^^(o.^y), i? a Gff 1 (o I r;ff 1 (r c )), uef^tjW), 

X G L 2 (0, T; H 2 {T C )) D L°°(0, T; ff 1 ^)) n fl^O, T; L 2 (r c )) 
for all T > 0, fulfilling the initial conditions 

tf(0) = tf° a.e. inn, i9 s (0)=i?° £ a.e. inr c (4.18) 
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and (|2.25p - (|2.26p . and the equations 

e \ $ t v+ / <9tln M (#)<u - / div(u t )v + e / Vi?t Vt> + / WVt> 

+ / fc(X)(0-0 s )t; = (M) VveV a.e. in (0, +00) , 

e! d t $ s v+ I d t ln^ s )v + e ! Vd t # s Vv 
jy c «^r c j r c 

- / d t \(X)v+ j W s Vu = j k(X)(0-0 a )v VwGff 1 ^) a.e. in (0, +00), 

6(u t ,v) +o(u,v) + / tfdiv(v)+ / (p M (X)u + 0^(11)) ■ v 

= w ,(F,v) w VvGW a.e. in (0, +00) , 
X t -AX + /3 M (X) + a'(X) = -A'(X)0 a - ^|u| 2 H M (X) a.e. in T c x (0, +00) , 

9 ns ^ = a.e. in 5r c x (0, +00) . 

Remark 4.3. Note that the approximate system (|4,19|) - (|4.22|) presents fewer technical difficul- 
ties than the analogous approximate version introduced in 0]. This is mainly due to the fact 
that, as we deal with the specific choice of ln(i?) and ln(i? s ) in the entropy (cf. Remark I Lip , we 
are allowed to directly introduce the Yosida regularization of the logarithm, instead of the more 
intricate approximating procedure exploited in [3]. 

4.2.1 Outline of the proof of global well-posedness for Problem (Pe). 

Since Problem (Pe) only slightly differs from the approximate problem considered in [4j, the 
global well-posedness for Problem (Pe), on any interval (0, T), can be obtained arguing in the 
very same way as in [U Sec. 3], to which we refer the reader for all details. Here, we shall just 
sketch the main steps of the proof. 

Step 1. First of all, one proves the existence of a local (in time) solution to (the Cauchy prob- 
lem for) Problem (Pe ) with the use of a Schauder fixed point argument. This involves 
establishing intermediate well-posedness results for each of the approximate equations, 
for which we refer to the calculations developed in [JJ Sec. 3.2], and defining a solution 
operator which complies with the conditions of the Schauder fixed point theorem (cf. [U 
Sec. 3.3]). 

Step 2. Next, to extend the local solution to the whole interval (0, T), one needs global (in time) 
a priori estimates. The latter substantially coincide with the ones formally performed in 
the proof of Proposition 12.51 and shall be repeated on the approximate system (|4,19p - (l4.22p 
within the proof of Proposition 14.51 

Step 3. Finally, uniqueness of solutions to Problem (Pe ) follows from the very same contraction 
estimates performed in [U Sec. 3.5]. 

Remark 4.4. We briefly justify our construction of the approximate Problem (Pe), referring 
to [U Rems. 3.2, 4.2] for more comments. 

The viscosity terms eJl^t) and e^Ri (dt$s) have been inserted in (|4.19p and (|4.2U|) . re- 
spectively, for technical reasons, related to the fixed point construction of a local solution for 



(4.19) 

(4.20) 

(4.21) 
(4.22) 
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Problem (Pg). In particular, the contributions eOlO&t) and e"Jl\ r (dt'&s) are essential to prove 
uniquen* of solutions to (the Cauchy problems for) approximate equations (|4.19p and (|4,20p . 
Furthermore, they also play a crucial role to make the estimate leading to the further regu- 
larity (|2T33aD - (12.43bl) rigorous, see the ensuing proof of Proposition 14.51 For the same reason, 
we have replaced the maximal monotone operators In, a and j3 by their Yosida regularizations, 
and correspondingly substituted the coupling terms X + u and — l/2£|u| 2 in equations (12.291) 
and (|2.31|) . with p fl (X)u and — l/2H M (X)|u| 2 , respectively. 

Now, as in the approximation of the system considered in [3], we shall keep the viscosity pa- 
rameter e distinct from the Yosida parameter p, in both approximate equations (|4.19p and (|4,20p . 
Thus, we shall prove the existence of solutions to Problem (P) by passing to the limit in Prob- 
lem (Pg) first as e \ for p > fixed, and then as p \ 0. This procedure shall enable us to re- 
cover on any interval (0, T) the L°°(0, T; H) -regularity for ln(i?) (the L°°(0, T; L 2 (r c ))-regularity 
for ln(i? s ), respectively) by testing the ^-approximation of (|2.27p (of (|2.28p . respectively), by 
the term ln^t?) (ln M (t? s ), resp.), and obtaining some bound independent of the approximation 
parameter p. In fact, such an estimate may be performed on equation (I4.19P (on (I4.20p . resp.) 
only when e = 0. For, if one kept e > 0, one would not obtain estimates on ln^(i?) independent 
of the parameters e and p, essentially because the term (e3£(i? t ), ln^(t?)) ( (e3£| r (dt'&s), 1% (i? s )) , 
resp.) cannot be dealt with by monotonicity arguments. 

4.3 Estimates on solutions to Problem (P^) 

Proposition 4.5. Assume (l2~HTl) - (l2~H5|) and (I2~H9]) - (I2.H12I) . Then, 

1. there exists a constant Kq > 0, only depending on the functions X, k, a, and on the 
quantity M (I2.4ip . such that for all e, p > the following estimates hold for the family of 
solutions {(t? e ^, "!?s j£ ^, u e ^, X £fl )} £tfl to Problem (P£); 



2. Furthermore, for all 5 > there exist constants Kf(5), K%(5,p) > 0, depending on 5, on 
the functions X, k, a, and on the quantity M (|2.4ip (K%(5,p) on p € (0,2) as well), but 
independent of e > and p > 0, such that the following estimates hold 



l|Vl?£ At ||L2(0,+oo;H) + ll V $s,£/JlL 2 (0,+oo;L 2 (r c )) < K 6 , 
e 1//2 ||^eMlli°°(0,+oo;y) + el//2 ||^s,£Mlli oo (0,+oo;H 1 (r c )) - K ® > 
W&enW L 00 (p,+oo;L 1 (Q)) + ll^*,6Mlli ao (0,+oo;i 1 (r c )) - K 6 ' 
||#eju - ^s,£Mlli 2 (0,+oo;L 2 (r c )) < ^6 , 

ll^t u eMlli 2 (0,+oo;W) + ll U £At (0,+oo;W) + ( U £At) IU°° (0,+oo) < ^6 
ll^eAtlli^CO.+oojH 1 ^)) + ll'9t ;><: £/i||L 2 (0,+oo;L 2 (r c )) 

+ ||/3/i(X £M )||2/x>(0 i+oo; £l(r c )) + ll H At (X eAt )||L oo (r c x(0,+oo)) < ^6, 
||9t ln M (l? e ,i)||i2 (o,+oo;V') + ll^t l n /i (^s,^) ||L2( 0)+oo; .ffi(r c )') < K 6 . 



(4.23a) 
(4.23b) 
(4.23c) 
(4.23d) 
(4.23e) 



(4.23f) 



(4.23g) 



\Wefi\\L°°(8,+oo;V) + ||$a,e/4 (5,+oo;_H' 1 (T c )) 

+ l|X e /J|z,°°(,S,+oo;.ff 2 (rc)) + ll^^e/i Hi 2 ((S.+oojH^r^JnL 00 (<5,+oo;L 2 (r c )) 
+ ||/3At(' )(: £A')lli <x> (5,+oo;i 2 (r c )) + ll u e/i llw 1 '° (<J,+°°;W) ^ ^(^) , 

|<9tt? eAt |[ L 2( 5>+00;i i2/7( n )) + ||dt#s,e/J|z, 2 («5,+oo;L 2 -p(r c )) < K 8 (5,p) , 



(4.24b) 



(4.24a) 
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Proof. We shall prove (|4.23p and (|4.24p by performing on Problem (P £ ) the same a priori 
estimates as in the proof of Proposition 12.51 and obtaining bounds independent of e, [i > and 
of t G (0,+oo). 



First estimate. We test (OS)) by tf £(U , (jOUl) by dOTI) by d t u £fl and by dtX £/J , add 

the resulting relations and integrate them on the interval (0, t), with t E (0, +oo). Basically, we 
put forth the same calculations as throughout (|3.1|) - (|3.7|) . up to the following changes. Instead 
of the formal identities (|3.ip . using the definition (j4.4j) of Jf^ we find 

o Jn Jo Jn 

■z^OM*))- I 

>Ci||i? eM (t)|| L i ( n)-C 2 -2||i?o||Li(n) ) 
the last inequality ensuing from (f47b|) and (|4.13bp . In the same way, we have 

/ / dtln^^)^ > Ci||0., e/1 (t)|Ui C r <O -C 2 - 2||^|ki ( r c) , 
Jo Jr c 

whereas we estimate 

in view of (|4.13a|) . and analogously for $ Si£ u- Then, we perform the same passages as in (|3.2|) — 
(13. 3D . up to replacing 3, (-) + , and /? with their approximations S M , p M , and Repeating 
the very same calculations as in (|3.6p - (|3.7p and taking into account some cancellations as well 
as (j4.13aj) and fl4.14aj> . we arrive at 

+ IWeMWm + \ I W^sXh + f [ k(X E „)(^ - 

z Jo Jo Jv r 



2\\ u eti\ L )\\V 



>0 JO JT C 

,eM\m{T c ) + \\Mt)\\LHF c ) + f l|W s " 2 
ft 



2 II ^^mv /ii-ff 1 (r c ) ' ll sV / Mi 1 (re) ' / II v w s\\L^(r c ) 

J 

+ C b f^ \\dt^A + ^K„(*)||w+ a^(t)) + lf r Pm(M*))IM*)| 2 

+ ^ ll^ll^(r c ) + ^l|VX^(t)||i 2(rc) + ^ (^(X £M (t))+a(X £M (i))) 

< c(l + ||0 o ||Li(n) + \\MW> + lKlki(r c ) + ll^lki(r c y 
+ ll u o|lw + «( u o) + l|VXo|lz,2(r c ) 

+ ll x o||L2(r c )l|uo||w + \W( x o)\\mr c ) + \W(^o)\\mr c ) + I \\ h \\v> 

Jo 

+ ll F lll°°(0,+oo;W')) + iQil/2 / H^MI^Hi 1 ^) + / ll F illw'|l u £Mllw • 



As in the proof of Proposition [231 from this inequality we deduce the bounds (|4.23ap - (|4.23c|) and 
(I4.23cl) - (l4.23fp (note that estimate (13. 8|) holds on this approximate level, too, thanks to (|4.9|) . 
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and the estimate for H fl (X £fl ) trivially follows from (|4.1ip ). Further, (|4.23dj) ensues from the 
bound 

||(fc(M) 1/2 (^M " <W)IU 2 (0,+oo;L3(r c) ) < C (4.25) 
and from ()2.H5j) . Again, the bound Q4.23gP for ln M (?? £/i ) and dt ln M (i9 S!£At ) is a consequence of 



the previous estimates and a comparison in equations (|4.19p and (|4.20p . 

Second estimate. We test KT9\i by tanh(-)dt$ £At , <^20\\ by tanh(-)dt# s , £M , differentiate (l43Tj) 
w.r.t. time and test it by tanh(-)9tu e/i , and differentiate (j4.22j) and multiply it by tanh(-)^X e/J . 
Notice that such an estimate is rigorous on this level, for dt~& £ ^ £ V and <9fi9 Sj£M € if 1 (r c ). 
Further, we sum the resulting relations and integrate in time. In place of using the formal 
identities (|3.13p . we have the following inequalities, which are direct consequences of (I4.6bj) : 

d t ln (U (i? £(U (-, r)) tanh(r)5 t i? eAt (-, r)dr 



o Jn 

t 



tanh(r)|c^ £M (-, r)\* In' (^(-, r))dr 

o JQ 

/ / tanh(r) ,'^^' r)|2 dr 



(4.26) 



4 T [ tanh(r)|,9 t GV 2 (-,r)| 2 dr, 
Jo Jn 



'o Jn 
and, analogously, 

/"* / ^ln M (i? w (-,r))tanh(r)^ s , eM (-,r)dr> /* / tanh(r)|^G]/ 2 (•, r)| 2 dr , (4.27) 
io ir c Jo ir c 

where we have set for all e, [i > 

e £M := |^ £M |+/i + 2, G s , £/X := |tf s , £M | + ^ + 2 . (4.28) 

Then, we observe that the remaining computations are not affected by the factor tanh(-) in a 
substantial way. In fact, the same calculations as in the proof of the second (formal) estimate 
in Proposition 12.51 go through, up to dealing with the integration by parts more carefully due to 
the presence of tanh(-). In particular, we have 

I [ W £M (-,r)V(tanh(r)<9ttf £M (-,r))dr 
Jo Jn 

\ f tanh(i)|V^ £M (t)| 2 -i f tanh'(r) f | W £/ ,(-, r)| 2 dr (4.29) 
1 Jn 1 Jo Jn 



2 

^ tanh(i) 



/ n iv^(t)| 2 -^, 



in which we have used the fact that tanh(0) = 0, that < tanh'(r) < 1 for all r G R, and 
estimate (|4.23ap . We handle the corresponding term for $ Sj£/i in the same way. Next, we 
estimate 

r fc(X CM (-, r)) tanh(r)(^(-, r) - (^(-, r) - $ s ^{; r)) dr 



o Jr 



(4.30) 



tanh(t) 



k(X e „(t)) (0 eit (t) - ^, £At (t)) 2 + ho + hi 
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where, integrating by parts 



1 



ho = -7 i tanh'(r)k(X £fl (;r))($ £ll (;r)-$ s ^(;r)y>-C (4.31) 

1 Jo Jr c 

J n = -\l I tanh(r)fc'(X £/U (-, r))d t X £fM (;r) r) - s>e|t (., r)) 2 dr 



2 . 
~ 2~ 



o Jr 



o 



I^X £A1 || L 4 (rc) ||i? £(U - $ s , £At ||i,4(r c ) ||# e/ * - ^,^llL 2 (r c ) ( 4 - 32 ) 



the inequality in (|4.3ip due to (I4.25|) . while the first passage in (|4.32p ensues from (l2.H5h and 
the second one from the Holder and Young inequalities, v being a positive constant to be chosen 
small enough. Furthermore, also taking into account (j2.5[) . we have 

tanh(r)6 (<9 t 2 u £(U (-, r), d t u £fJ ,(-, r)) dr 

> ^tanh(t)||^u £M (t)||^ - ±j\ a nh>(r)b(d t u £ll (;r),d t u £ll (;r))dr (4-33) 

>^tanh(t)||^u £M (t) 
the last passage due to (|4.23e|) and (|2.5p . In the same way, we find 

j\ a nHr)dlX £ ,(-,r)d t X £ ,(;r)dr > t -^l\\d t X £ ,(t)\\ 2 LHrc) - *± , (4.34) 
tanh(r)V(5 t X £M (-,r))V(a^(-,r))dr= /"\anh(r) / |d t (VX £M (-, r))| 2 dr , (4.35) 



'o Jv c 
and we remark that 

■■I 



tanh(r)(a' (u e J-,r))d t u E J-,r),d t u E J-,r))dr > 0, 



o 



(4.36) 



tanh(r) / #,(X B/1 (-, r))|d t X £/i (-, r)| 2 dr > . 



Finally, we perform the same computations as in (I3.19p - (l3,20p . and observe that 

I I tanh(r)a t ( PjU (X £/i (-,r))u £M (-,r))9 t u £M (-,r)dr = /i2 + /i3, (4.37) 
jo Jv c 

h2 = [ I tanh(r)p M (X £M (^r))|d t u £/ ,(-,r)| 2 dr >0, (4.38) 
J o Jr c 

^13= / tanh(r) f H M (X £/i (-, r))<9jX £M (-, r)u £fl (-, r)d t u £tl (; r)dr 
Jo Jr c 

[* 3 
> — C / ||9tX £ ^||x,2( rc )||u £/L4 ||w||^U £/i ||w > — C^6 j 
JO 



(4.39) 



the passages in (|4.39p due to the Holder inequality, the continuous embedding (12. 2p . and esti- 
mates (|4.23ep - (j423f|) . In the end, we point out that the analogues of (pnU]) - (PTTT|) and (pTT5|) - 
(13.161) go through. Then, collecting (|4.26p - (l4.39p and putting forth the same arguments as in 
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the proof of Proposition 12.51 we conclude estimate (|4.24ap for all 5 > 0. We also find 

\\dt@l(?\\L 2 (6,+oo;H) + |l^t©s^lU 2 (<5,+oo;i 2 (r c )) < C<5 > 

and, arguing in the same way as in (|3.24p - (|3.25p . we infer that for all p G (0, 2) 

|l^t@£Ai||L2(5,+oo;i 12 /7(Q)) + ll^@s,£Mlli 2 (<5,+oo;L 2 -P(r c )) < C S,p , 

whence the bound (|4.24bp ensues. n 
4.4 Proof of Theorem [231 

4.4.1 Passage to the limit in Problem (Pe) as e \ 

First, we introduce the boundary value problem obtained by taking e = in Problem (Pe), 
which we shall supplement with the initial data (|2,lip - (|2.14p . 

Problem (P M ). Given a quadruple of initial data (i?o> $s > u 0> <^o) complying with ()2.11 j) (j2. 141) . 
find functions (i?,i? s , u,X), u with the regularity (j2TTHD . X with (j23U1) . and, for all T > 0, 

G L 2 (0,T;y) nL oo (0,T;L 1 (^)), ln^(^) € -ff 1 (0, T; V') , (4.40) 
^ s eL 2 (0,T; J ff 1 (r c ))nL oo (0,T;L 1 (r c )), l nfl (^)eff 1 (0,T;ff 1 (r c )') 1 (4.41) 



(4.42) 



fulfilling the initial conditions 

ln M (0(O)) = M^o) a.e. inO, M^( )) = In M (^) a.e. inr c , 
u(0) = uo a.e. in 0, X(0) = Xq a.e. inT c , 

the equations 

(d t \nJ#),v) - [ div(u t )v+ / WV« + [ k(X)(&-0,)v 

Jn Jn Jr c (4.43) 

= (h, v ) Vw € V a.e. in (0, +oo) , 

{d t ln^ s ),v) - I d t \(X)v+ [ V# s Vv 
Jv c Jr c 

= [ k(X)(0-0 a )v VoGF 1 ^) a.e. in (0, +oo) , 
Jr c 

and such that (u,X) comply with (|43I^ - (|Q5|l . 

Notice that no uniqueness result is available for (the Cauchy problem for) Problem (P M ). 
The following result is a straightforward consequence of Proposition 14.51 

Proposition 4.6. Assume (f2~HT|) - pdTo]) and (j2Jl9]) - (|2.H12p . Let /i > be fixed. Then, there 
exists a (not relabeled) subsequence of {(t) Efl ,'& S:e ^,u Efl ,X £ ^)} E and a quadruple ($,'!? s , u, X) such 



(4.44) 
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that for all T > the following convergences hold as e \ 

ti £ll -n? in L 2 (0, T; V), e%(-& eil ) ^ in L°°(0, T; V 7 ), (4.45a) 

0,^^-8, inL 2 {Q,T;H 1 {Y c )), e#r c (tf s , £M ) - m L°°(0,T; iJ 1 ^)'), (4.45b) 

h V (tf £/ J - ln^tf) mL 2 (0,T;7), (4.45c) 

ln M (^ £A1 ) -> ln M (i9) in C°([0, T];X) for every Banach space X with V'il, (4.45d) 

eK(i? eM ) + WW *« H H^ T; V) , (4.45e) 

M'W M*-) w ^ T ; ffl ( r c)) , (4.45f) 
liV(^, £M ) -» ln M (i?«) in C°([0, T]; F) /or euery 5anac/i space F mtt iT 1 ^)' <s F , (4.45g) 

^(W + HtW^M^) inff 1 (0 j r;H 1 (r c ) / ), (4.45h) 

X £M -*X in L 2 (0, T; # 2 (r c )) n L°°(0, T; -fT 1 (r c )) n H\0, T; L 2 (r c )), 

X £At ^X ^nL 2 (0,T;F 2 -^(^ c ))nC ([0,T];F 1 ^(^ c )) for all p G (0, 1), 

p M (X £M )-*p M (X) zn L°°(0, T; tf 1 ^)) n fl^O, T; L 2 (r c )), 

p M (X £M ) -> p M (X) mC^O.rljtf-^rj) /or p G (0, 1), 

/MM^/W mL 2 (0,T;L 2 (r c )), (4.45k) 

H M (X £(U )^*H M (X) mL°°(0,T;l7(r c )) /or a// 1 < p < oo, (4.451) 

a M (u £M ) a„(u) in L 2 (0, T; W), (4.45m) 

u £M ^u mff^TjW) 

u £M ^u mO^O^iT 1 ""^) 3 ) for all p € (0, 1) , 



(4.45i) 
(4.45j) 



(4.45n) 



(the symbol <s in (]4.45d|) and ( 4.45gP signifying compact inclusion), and (i?, i? s , u, X) is a solution 
to Problem (P^). 

Furthermore, (i?, i? s ,u, X) /ios i/ie additional regularity (|2.36aj) (|2.36cj) and (|2.36e|) on e?;e- 
ry interval (6,T), for all < 5 < T and, up to the extraction of another subsequence, for 
{(i? £A j, "!? Sj£At , u £M , X £/J )} £ the enhanced convergences hold as e \ for all < 5 < T 



(4.46a) 

(4.46b) 
(4.46c) 
(4.46d) 



tf £M -* in L°°(5, T; V) n H l (5, T; L 12 / 7 (0)), 
^ ^ # in c°([5, T]; H 1 ~ P (Q)) for all p G (0, 1), 

inL°°(,5,r;F 1 (r c ))n J ff 1 (^r;^ 2 ~ p (r c )) for all P e (0,2) 

$s,en ^ in C°([5,T]; iy 1 "'^)) /or a// p G (0, 1), 

u £/ ^*u mW^ftTjW), (4.46c) 
X £M -* X zn L°°(5, T; tf 2 (r c )) n H\5, T; H l (T c )) D ^• 00 (5, T; L 2 (r c )), 
X £ ^X znC°([5,r];F 2 -"(r c )) for all p e (0,2), 

/3 At (X £/x )-*/3 /i (X) inL°°(,5,r;L 2 (r c )). (4.46e) 

Sketch of the proof. Since the proof of this result substantially goes along the same lines as 
the proof of [U Prop. 4.4], we shall just outline its main steps, referring to [I] for all details. 

• First of all, using estimates (|4.23p and standard compactness results (cf. |18j). with a 
diagonalization process we extract a subsequence of {($ £/i , i? s ,e/i) u e^' ^e^)}e for which con- 
vergences (|4.45ap - (|4.45np hold as e \ 0. In order to show that •& G L°°(0, T; L 1 ^)), as 
in the proof of [U Thm. 1] we exploit a Lebesgue point argument. Indeed, using the first 



Bonetti, Bonfanti, Rossi / Analysis of a thermomechanical model for adhesive contact 33 



of (|4.45ap and the convexity of | ■ |, we find that for all to € (0, T) and r > such that 
(t -r,t + r) C (0,T) 

/ ° + / \0\ < liminf ( [ ° + [ \0 elt \) < 2rsup ||0 B/ ,|U<»(o,T;Li(n)) < 2rK 6 , (4.47) 
Jto-r Jn £ \° \Jto-r JQ. J e>0 

the latter inequality due to (|4.23c|> . We now divide the above relation by r and let r J, 0. 
Using that the Lebesgue point property holds at almost every to G (0, T), we obtain for 
all T > the estimate 

l|0|U«(O,T;£i(n)) < *6 ■ (4-48) 
A completely analogous argument may be developed for proving that $ s 6 L°°(0, T; L 1 (r c )). 

• The identification of the limits of the sequences {/^ (X e/J )} and {H M (X £M )} follows from the 
strong- weak closedness of the graphs of the operators /3u, and H^. As for the limits of 
{a^(u e/ J}, {ln M (i? £/i )} and {ln^^^)}, thanks to [U Lemma 1.3, p. 42], it is sufficient to 
prove that for all t > 

limsup T / lrn^)^ < f / ln M (tf)tf, 
e\o Jo Jn Jo Jn 

limsup / / ln M (^ S)£/1 )^ S)£jU < / / ln^i?*,)^ , 
e\o Jo Jr c Jo Jn 

limsup / / a^Uep) • u £(U < / / o^(u)-u. 

e\o Jo ir c jo Jr c 

Convergences ( |4.45j[ ) for {p M (X £M )} ensue from the corresponding (|4.45ip for {X £ ^} and 
fron the fact that < p^(r) = H M (r) < 1 for all r G R, cf. (jlTD) . Notice that the pair 
(u, X) complies with the initial conditions (|2.25p - (|2.26p in view of convergences (|4.45np 
and (|4.45ip . In the same way, combining convergences (|4,13dp and ()4.14dp of the sequences 
{ln/t(#e)}> { ir v(^s e )}> with convergences (I4.45dp and ( |4.45gP , respectively, we conclude 
the initial conditions (|4.42p for ln^(i?) and m^(i9 s ). Finally, we refer to the proof of [H 
Prop. 4.4] for the argument showing that $ s ,u, X) is a solution to Problem (P^). 

• The further regularity ()2.36ap ~ (|2.36cl) and (I2.36ep for (■&, i9 s , u, X) and the enhanced con- 
vergences (|4.46p ensue from estimates (|4.24p , see also the proof of Theorem [ 



n 

In view of the definition of approximable solution to Problem (P) which we shall give in Sec- 
tion we select, among solutions of Problem (P^), only the ones just constructed by passing 
to the limit in Problem (P £ ) as the viscosity parameter e vanishes. Albeit improperly, within 
the scope of this section we shall refer to them as viscosity solutions. 

Definition 4.7 (Viscosity solutions of Problem (P^)). Let fi > be fixed. We say that a 
quadruple (i?, i9 s , u, X) is a viscosity solution of Problem (P M ) if 

1. it is a solution to Problem (P^); 

2. there exist a sequence Ek \ and a family of solutions {(i9 £fe/J , i? Sj£fe/J , u £fe/J , X £feAt )}^ of 
Problem (P £fe ) converging as \ to i? s ,u, X) in the sense specified by (|4.45a|) - 
(I4.45np . on every interval (0, T). 

In particular, it follows from the last part of Proposition 14.61 that every viscosity solution to 
Problem (P^) has the further regularity (I2.36ap - (|2.36c|) and (|2.36e|) . 
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4.4.2 Passage to the limit as p \ 

Proposition 4.8 (Estimates on viscosity solutions of Problem (P^)). Assume (|2.H1|) (|2.H5|) 
and (pEH9l) - (l2.H12p . Then, 

1. estimates (|4.23|) on the half-line (0, +oo), (|4.24a|) on (5, +oo) for all 6 > 0, and (|4.24bp 
on (5, +oo), for all 5 > and p G (0,2), hold, with the same constants Kq, Kf(8), and 
K$(5, p), for all p > and for every viscosity solution ^ s ,^, u M , X^)} to Problem (P^). 

2. Furthermore, for all T > there exists a constant Kg(T), only depending on T, on the 
quantity M (I2.41|) . on the functions A, k, anda, as well as on || ln(i?o)||,tf and || ln(i9g)||£2(r c ), 
such that for all p > 

II m /i (^)|| L o O(0; r;H) + II M0 a>M )|U»(o,Tj£3(r c )) ^ Kg ( T ) • ( 4 ' 49 ) 

As already mentioned in Remark 14.41 the calculations leading to (I4.49|) cannot be performed on 
the solutions of Problem (Pe ), with e > 0. 

Sketch of the proof. The first part of the statement is an obvious consequence of the Defini- 
tion (14. 7p of viscosity solution, of estimates (|4.23p and (I4.24p . and of a trivial lower-semicontinuity 
argument. The enhanced estimate (I4.49P for in. fl ('&u,) and ln^(-d SjfJj ) is obtained by testing (14.431) 
by ln^i?^), (|4.44|) by ln^i?^), adding the resulting relations and integrating on time. We refer 
to [U Sec. 4.2] for all details. Here, we just point out that, being ln($o) £ H, thanks to H4.42H 
there holds 

|| hvOV(0))||tf = || ln^ )\\ H < || ln(tf )||tf , 

and the analogous bound for || ln /Lt (i? SjAi (0)) ||i2(r c ) ensues from the fact that ln(i?°) € L 2 (F C ). 
Such estimates are then used in the calculations yielding (I4.49p . □ 

Thus, we have the following proposition, which concludes the proof of Theorem 12.11 

Proposition 4.9. Assume ([2TH11 - (|2?H5l) and pH9|l - (|2.H12jl . Let {(^, # SjM , u M , X^)}^ be a 
sequence of viscosity solutions to Problem (P M ). Then, there exists a (not relabeled) subsequence 
°f ^s,a, u fj,i ^-u,)}^ and functions ('&,'& s ,u,X,r),£,Q such that for all T > the following 
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convergences hold as n \ 

in L 2 (0,T;V), (4.50a) 

&s,v^&s mi 2 (0,T;F 1 (r c )), (4.50b) 

liV(^H*ln(??) in L°°(0,T;H)nH 1 (0,T;V'), (4.50c) 

ln M (# M ) -> ln(#) in C°([0, T]; X) /or ewer?/ Banach space X with H m X , (4.50d) 

ln M (^ jM )^*ln(^) inL^TsL^r^nid^T^rc)'), (4.50e) 
ln M (tf SjM ) -> ln(0 s ) in C°([0, T];Y) for every Banach space Y with L 2 (T C ) m Y , (4.50f) 
X^*X tn L 2 (0, T; H 2 (T C )) n L°°(0, T; tf 1 ^)) D H^O, T; L 2 (r c )), 
X^X znL 2 (0,T; J ff 2 ^(r c ))nC ([0,T]; J ff 1 ^(r c )) /or aHp G (0,1), 
P/AH** + in L°°(0, T; ff 1 ^)) n ^(0, T; L 2 (T C )), 
M X »)^ X+ tnCV^T^H^irc)) for all p G (0,1), 

/SmW- 1 ^ inL 2 (0,T;L 2 (r c )), £ G /3(X) a.e. mT c x (0, T), (4.50i) 
H M (X M )^*C znL°°(0,r;LP(r c )) for all 1 < p < oc, 
C G H(X) a.e. inT c x (0,T), 

a M (u M )^r/ in L 2 (0,T;W), rj G a(u) a.e. in(0, T), (4.50k) 
u^u in fl' 1 (0,T;W) 

u^u zn C ([0,T];i7 1 -^(O) 3 ) /oraHpG(0,l), 



(4.50g) 

(4.50h) 
(4.50i) 
(4.50j) 
(4.50k) 
(4.501) 



and i? s , u, X, T7,£, £) is a solution to Problem (P). 

Furthermore, the functions #, i? s , u, X, and £ /iawe t/ie additional regularity (12.360 on e?;e- 
ry interval (5,T), for all < 5 < T and, up to the extraction of another subsequence, for 
$s,/ij u ^t) X^)}^ ^ e enhanced convergences (I4.46P /ioW as /i \ /or allO < 5 <T. 



Sketch of the proof. The proof follows the very same lines of the argument for Proposition 
(cf. also the proof of [U Thm. 1]). We just point out that convergences (|4.50hp ensue from the 
pointwise convergence p At (X At ) — ► X + a.e. in T c x (0, T) (which can be deduced from ( |4.50gP 



exploiting formula (j4.10|) ). and again from the bound (|4.1ip on p' . Furthermore, (|4.50cp ~ (l4~50fl) 
are a consequence of the additional estimate (|4.49p . Finally, as in the proof of Proposition l4.6l the 
strong convergences (I4.50d|) and (14.50f|) . as well as the Yosida convergences ln^(i?o) — l n (^o) 
in H and ln^(i?^) — 1 ln(i?g) in L 2 (r c ) as [i \ 0, enable us to pass to the limit in initial 
conditions H4.42H and deduce that 



C ! 



ln(tf(0)) = ln(tfo) a.e. in O, ln(tf s (0)) = ln(t?°) a.e. in T c 
whence initial conditions (|2.23p and (|2.24p . □ 

Remark 4.10. Notice that, in the statement of Proposition (14. 9p . assumptions (I2.H10|) - (I2,H12|) 
on the data f, g, and h are stronger than the data requirements in Theorem 12. 1\ and there is the 
additional condition (|2.H9p on (3. This is due to the fact that, to avoid unnecessary repetitions, 
we have chosen to unify in the present section the proof of global existence with the proof of the 
long-time estimates. 

In fact, a closer perusal of the proof of Theorem 12.11 and a comparison with the argument 
for |H Thm. 1] show that the sole ([2~HT]) on and (f2~H6l - ([2~H8|) on the data h, f , and g are 



Bonetti, Bonfanti, Rossi / Analysis of a thermomechanical model for adhesive contact 36 



sufficient for the global existence of solutions to Problem (P) on the finite-time interval (0, T). 
Further, under assumptions (|2.35p (which are the finite-time versions of (|2.H10p - ([2.H12p ) one 
proves the further regularity (|2.36|) on the interval (0, T). In particular, condition (|2.H9j) on (3 
is not necessary for the aforementioned finite-time results. 

4.5 Rigorous proof of Proposition 12.51 

We are now in the position of specifying approximable solutions to Problem (P) as the ones 
arising as limits of viscosity solutions to Problem (P^). 

Definition 4.11 (Approximable solutions of Problem (P)). We say that u, X, r], £, £) is 

an approximable solution of Problem (P) if 

1. it is a solution to Problem (P), 

2. there exist a sequence \ and a family of viscosity solutions {(i?^, "&s,ii k i u /j, k i ^W)}* 
of Problem (P Mfc ) converging, as \ 0, to (i?,i? S) u, X,rj,^,Q, in the sense specified 
by (|4.50ap - (|4.501|) . on every interval (0, T). 

Remark 4.12. Taking into account Remark 14.101 we might state Theorem 12. II in the following 
more precise way: under assumptions (|2.Hip - ([2.H8p . on every interval (0, T) Problem (P) admits 
at least an approximable solution $ s , u, X, rj, £, £). If, in addition, (|2.35|) holds, then every 
approximable solution has the further regularity (|2.36|) . 

Rigorous proof of Proposition 12.51 It follows from the Definition (|4.11|) of approximable 
solution, from the estimates on viscosity solutions to Problem (P^) specified in Proposition 14.81 
and from elementary lower-semicontinuity arguments. n 
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